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carried out to complete the description of jet modification at high energies at future high energy 
colliders. 
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1 Introduction 



During the past nine years, high energy heavy-ion colhsions have been studied at the the Relativistic 
Heavy-Ion Colhder (RHIC) at Brookhaven National Laboratory (BNL). The three most striking findings 
in heavy-ion colhsions at RHIC are the observation of a large azimuthal asymmetry in soft [pt < 2 
GeV) hadron production, denoted as elliptic fiow; the observation of a large suppression of hadron yields 
at high transverse momentum [p^ > 6) GeV and the finding that the yields and the elliptic fiow at 
intermediate pt follow a scaling behaviour that separates baryons and mesons by a factor of 3/2. Each 
of these observations indicates that hot and dense strongly interacting matter is formed in collisions at 
RHIC, possibly a deconfined Quark-Gluon Plasma. Soon, the LHC will start colliding nuclei at much 
larger energies, which will provide critical tests of the current understanding of hot and dense QCD 
matter at RHIC. 

In this review, we will concentrate solely on the suppression of high-p^- hadron production and 
specifically its theoretical description within a perturbative QCD (pQCD) based formalism. In the 
collision of two heavy-ions at very high energies, most of the hard valence partons go through the 
collision un-defiected. The prevalent interaction is that between the softer "sea" partons. It is these 
soft interactions which lead to the formation of the hot and dense matter. Here "soft" means involving 
momentum scales of the order of Aqcd- Occasionally the hard partons undergo hard scattering (at 
scales much larger than Aqcd) leading to the formation of two back-to-back hard partons with a large 
Pt- High-pT hadrons originate in the fragmentation of such high px partons after their escape from the 
medium. The hard partons lose energy through interactions with the hot and dense medium. High pt 
signatures are of special interest, because they are expected to be described by perturbative QCD due 
to the hard scale of the jets. At the LHC, it will be possible to reach px ~ 100 — 200 GeV, i.e. a factor 
of 10 increase in the hard scale of the jets compared to RHIC. Jet physics based on perturbative QCD 
is therefore expected to play a central role in the study of hot QCD matter at the LHC. 

While the development of the theoretical framework to describe parton energy loss was started 
a little less than 20 years ago jl], the basic observation that high-p^- hadrons are suppressed was 
only made with the first RHIC data [21 [3]. Since then there has been significant progress in our 
understanding of parton energy loss. On the experimental side, the variety of measurements has greatly 
increased, now extending to a variety of di-hadron and heavy fiavour measurements. The precision 
and momentum reach of these measurements is also much more extended compared to the first results 
from RHIC. Phenomenological calculations have become much more detailed, incorporating a variety of 
initial state effects such as shadowing, Cronin broadening [1] and more realistic models of the collision 
geometry based on hydrodynamical calculations constrained to describe the soft hadron spectrum. 
There are now four different theoretical schemes for energy loss calculations, which incorporate different 
physical assumptions regarding the scales involved and the microscopic structure of the medium. These 
approaches can be compared and contrasted to provide insight into the essential aspects of the dynamics 
of energy loss in heavy ion collisions. 

In this review, we aim to provide a survey of the current state of the theoretical understanding of 
parton energy loss and the experimental tests of this understanding. Our objective is not to provide an 
exhaustive compilation of all experimental results and theoretical calculations that have been performed, 
but rather to provide an overview of the current understanding and methods which can serve as a 
'reading guide' for the extensive literature on the subject. The theoretical part of the review introduces 
the pQCD based formalism of jet quenching as an extension of the standard factorized approach to 
hard processes in vacuum and then discusses the similarities and differences between the four different 
theoretical approaches to parton energy loss in some detail. In the experimental part of the review the 
emphasis is on a subset of observables which are, or can be, rigorously calculated within pQCD and 
focus on what has been learnt from these. This field is under active development, both theoretically 
and experimentally. On the theory side, the topic of most current interest is the development of Monte 
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Carlo routines to address some of the shortcomings of the analytical approaches and to compute more 
exclusive jet based observables. On the experimental side, measurements have been extended to include 
photon-jet events and multi-hadron observables, including first attempts to reconstruct jets in heavy 
ion collisions at RHIC. In conclusion, we will discuss these developments and how they are expected to 
address some of the main open question about parton energy loss. 

In the remainder of this chapter we provide a very brief introduction of how the factorized formalism 
arises in pQCD. In Chapter 2, we re-derive the multiple scattering induced single gluon emission cross 
section within this language. This is most straightforwardly achieved in a particular variant of what has 
come to be denoted as the higher twist approach. After this we relate our results to the derivations in 
other schemes. In Chapter 3, we describe the inclusion of multiple emissions. This is as yet a theoreti- 
cally unsettled regime and we review the three different means by which multiple emissions have been 
included. This is followed by a direct comparison of the various formalisms, where the medium modified 
fragmentation function is calculated in an identical medium. In Chapter 4 we extend our formalism to 
include heavy flavor modification along with elastic loss and diffusion. In Chapter 5 we describe the 
phenomenological setup used to compute jet modification in a static large nucleus and in a dynami- 
cally evolving quark gluon plasma. In Chapter 6, we compare the results of recent calculations in the 
literature on single hadron, dihadron, photon triggered and heavy flavor observables to measurements 
in Deep-Inelastic Scattering (DIS) from HERA and those in high-energy nuclear collisions from RHIC. 
In Chapter 7, we provide an outlook to the LHC and forthcoming theory calculations. 



1.1 Background 

Quantum chromodynamics (QCD) is now the accepted theory of strong interactions. In spite of the 
large diversity of hadronic states observed in nature, the Lagrangian of QCD is rather simple, involving 
only two types of fundamental fields: quarks and gluons, interacting via an SU {Nc = 3) gauge theory [5] 

Uf 

^QCd{x) = -\f^,A^)F^"'{^) + E^'(^) h'^M,, + m,] rA^)., (1) 

q=l 

where i,j run from 1 to Nc representing the colors of the fundamental quarks of flavor q and a runs 
from 1 to A'^^ — 1 representing the colors of the adjoint gluon gauge field. For the purposes of this review 
Uf will mostly be limited to 3, for the light quarks u, d, and s, which will be treated as massless. The 
heavy flavors charm and bottom will be discussed separately and will have non-zero mass terms. The 
top quark will not be included in this review (thus Uf < 5). The covariant derivative and the adjoint 
field strength have the usual definitions, 

= kjd'' - ^9tl^A^'' and F^^^ = d^A^ - d.A'^ + gr'^A^Al. (2) 

In the equation above, is the gluon four-vector potential with adjoint color a, t^j are the Gell-Mann 
matrices, f"'^'^ is the completely antisymmetric tensor and g = -y/lvraj is the strong interaction coupling 
constant. In this review, we will often refer to the fine structure constant also as the "coupling". 
Which is meant, will be denoted using the appropriate symbol. 

Like all renormalizable quantum field theories, the coupling depends on the renormalization scale 
fi, usually chosen as the relevant hard scale in the problem to minimize the effect of higher order 
contributions. Unlike the case of QED or weak interactions, QCD is asymptotically free, which implies 
that the coupling becomes weaker as the scale is raised [6l [7] . Solving the (one loop) renormalization 
group equation, and comparing with experimental data to obtain the initial condition, one obtains, 

= T ^ 7 (3) 
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The above equation implies that at scales far above Ivqcd ~ 200MeV, QCD should be weakly coupled 
and a perturbative expansion in terms of should become feasible. The applicability of this statement, 
of course, depends on the process in question. In the case of the total cross section in e^e~ annihilation, 
single hadron inclusive e^e~ annihilation, the total cross section in Deep Inelastic Scattering (DIS), 
single hadron inclusive DIS and the Drell-Yan effect, pQCD begins to become applicable beyond a 
^ 2 GeV^ [H [9l [To]. The quantity here refers to the relevant invariant mass that sets the 
hard scale in the problem. In all the three cases mentioned, this is the invariant mass squared of the 
intermediate photon. In the case of hadron-hadron scattering, one requires a specific hard momentum 
transfer process, such as the production of a high transverse momentum (high pt) particle. The hard 
interaction in this case is mediated by the strong force (meaning the higher order corrections are of a 
different type) and such processes seem to need a minimum Pt 2 GeV (Q^ ^ 4 GeV^) for pQCD to 
be applicable |TT]. In all cases mentioned above the applicability of pQCD refers to the agreement of 
calculations based on pQCD and experimental data. 

With the diminishing of momentum transfers, the strong coupling fine structure constant as{Q'^) 
continues to grow and confines all particles which carry color charge within composite color singlet 
hadrons. While one may still construct an effective Lagrangian using the underlying symmetries of 
QCD and carry out a perturbative analysis, such theories will not constitute any part of this review. 
All interactions with a momentum scale below a ~ 2 GeV^ will be considered as non-perturbative and 
will usually be contained in a non-perturbative distribution such as a fragmentation function, a parton 
distribution function or an in- medium transport coefficient. The separation of the non-perturbative 
part from the hard perturbative part of the calculation will be discussed in the next subsection. 



1.2 pQCD and factorization 

The ability to apply pQCD to describe a particular process should not be confused with the application 
of perturbative expansions in other gauge theories such as QED, the electro- weak theory or even in low 
energy effective theories of QCD such as chiral perturbation theory. In those theories, the incoming 
and outgoing asymptotic states in a scattering event carry the quantum numbers of the fields in the 
respective Lagrangian densities. While the QCD Lagrangian is cast in terms of quarks and gluons 
and by pQCD we do mean a diagrammatic expansion involving those fields, asymptotic states in strong 
interactions, due to confinement, are never quarks or gluons but composite hadrons. The ability to apply 
pQCD thus means simply the ability to isolate a section of the interaction which can be systematically 
described using a perturbative expansion in as involving quarks and gluons, from the remaining part 
of the process which is non-perturbative. 

In all the reactions mentioned above where pQCD is applicable, there exist sub-processes over a 
range of energy scales up to the hard scale Q^. Most of these cannot be described using pQCD and 
need to be separated from those sub-processes which involve hard scales. The technical machinery which 
demonstrates this separation order-by-order in the coupling constant is called "factorization" [121 [131 IS] • 
The result of factorization is usually stated as a theorem with corrections power suppressed at very large 
Q^. We illustrate an example of this theorem for the case of a high energy p-p collision leading to the 
formation of a high pt hadron. Assume the collision is in the center of mass frame with each proton 
carrying a momentum P{—P). The physical picture of this process is one where a hard parton in one 
of the incoming nucleons, carrying a light-cone momentum fraction Xa {Pa = XaP^) scatters of a hard 
parton in the other nucleon with light-cone momentum fraction Xi, [p^ = —Xi,P~) and produces two back- 
to-back partonic jets. Light-cone momenta are defined as P+ = (P° + P^) /V2 and P~ = (P° - P^) /V2. 
The hadronization products of one of these jets (c) will yield the high pt hadron h carrying a momentum 
fraction z of the original jet's momentum. The factorization theorem for the differential cross section 
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(in pt and rapidity y) states that the process can be expressed as [13], [15] , 

- ' dXadXbGaiXa, IJ,f)Gb{Xb, fJ^f) 4 + O 7^ , (4) 



dyd^PT J dt nz \Q 

where Ga{xa)[Gb{xb)] is the parton distribution function (PDF) to find a hard parton with a momentum 
fraction Xa {xb) in the incoming proton, da/dt is the hard partonic cross section (with a partonic 
Mandelstam variable t) and D^(z) is the fragmentation function (FF), the distribution of hadrons with 
a momentum fraction z produced in the hadronization of the outgoing hard parton. The second term 
on the r.h.s. indicates corrections to the factorization theorem that are suppressed by powers of the 
hard scale Q^. The quantity A represents the scale of soft processes in the collision. This contribution 
becomes negligible as — oo. Even in this asymptotic limit, the only term in the equation above 
that can be completely calculated within pQCD is the hard partonic cross section; the PDFs and the 
FF are non-perturbative objects which represent physics at softer scales. 

The primary utility of factorization is based on the lack of interference between the hard and the 
soft scale, i.e., the PDFs and the FF (as well as the partonic cross section) are all squares of amplitudes. 
Eq. (jl]) contains only a convolution in two parameters {xa,Xb) between these probabilities. To compute 
the hard cross section one proceeds order by order and calculates all the amplitudes that may contribute 
to the process of two incoming partons scattering to n outgoing partons [n = 2 for Leading Order (LO), 
n = 3 for Next-to-Leading Order (NLO)]. These are then summed and squared to obtain da/dt. The 
non-perturbative quantities G{x,fi) or D{z,fi) cannot be calculated from first principles, however they 
have well defined operator expressions on the light-cone, e.g., the PDF to find a quark with momentum 
fraction x, in a proton state \p) traveling with a large light-cone momentum p"*", is defined as, 

G{x) = I ^e-^"^-(p|^(y-)7+^(0)b). (5) 

Given the operator expression, one may calculate higher order corrections. These tend to have 
collinear divergences which are absorbed into a renormalization of the operator product (or non- 
perturbative expectation values in the case of the FF). Such a redefinition introduces a scale dependence 
of the expectation. Divergent contributions up to this scale, denoted as the factorization scale /i/, are 
absorbed into the definition of the PDF. If fif is large compared to Aqcd, the change of the PDF 
and the FF can be calculated by means of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation [HI HH HB] , 



dG{x, jj) Us dy 



aiog/i 27: J ^ y \y J 

In the equation above, P{y) is the splitting function which represents the probability for a hard quark 
to radiate a gluon and still retain a momentum fraction y. The homogeneous DGLAP equation, as 
expressed in Eq. ([6]) apphes only to the non-singlet quark distribution [Gns{x) = Gq{x) — Gq{x)]. 
The singlet distributions {Gs{x) = [Gq{x) + Gq{x)]/2) have couplings that mix the quark and gluon 
distributions. 

Equation (|6]) is a differential equation and thus requires an initial condition for its solution. This 
requires the experimental measurement of the non-perturbative distribution at one value of the scale. 
Herein lies the other advantage of factorization: Once the non-perturbative distributions are factorized 
from the hard cross section, they are independent of the process and become universal functions in the 
sense that they may be given an identical definition and measured in a completely different process. In 
the case of totally inclusive DIS, the factorization theorem yields. 
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The operator expression for the PDF as well as its evolution equation, derived in this case, is identical 
to that in Eq. (jlj). Thus, the PDF measured in DIS (^^, may be directly substituted in Eq. (jl]) to 
compute single particle inclusive cross section in p-p collisions. 

To summarize, the applicability of pQCD to vacuum processes, such as in Eqs. (I4f7|) . consists of 
the ability to calculate the hard partonic cross section order by order, as well as to compute the 
scale dependence of the non-perturbative distributions. While these non-perturbative distributions 
need to be measured at one scale in experiment, they have rigorous operator definitions which are 
identical in all processes where they can be factorized and thus are universal functions. Even though 
the discussion above was focused on the PDF, an almost identical factorization theorem and DGLAP 
evolution equation may be written down for the FF. Most of the remaining review will deal with the 
vacuum and in-medium modification of the FF. 

1.3 Hard jets in Semi-inclusive DIS on a nucleus and heavy- ion collisions 

The preceding subsection briefly described the factorized formalism of high energy pQCD as applied 
to hard processes. The underlying assumption of factorization has been rigorously proven in all the 
examples mentioned above. In this subsection, we outline the extension of this factorized formalism 
to the modification of hard jets in a dense medium. Unlike the case of hard processes in vacuum, 
the factorization theorems which are the underlying assumptions in these calculations, have not been 
proven to the same degree of rigor. So far there have only been a handful of attempts in this direction 
both in pQCD [19l [201 HI] ^^11 as in an effective-field theory approach [22]. These factorization 
theorems will not be discussed further; we will assume their applicability. In this review we will focus 
on the application of the factorized theory to jet propagation in dense extended media and the ensuing 
phenomenology. 

Imagine the single-inclusive DIS of a hard lepton on a proton in the Breit frame. This is illustrated 
in Fig. [U The virtual photon with large negative light-cone momentum q = [Q'^/2q~, q~ , 0, 0] strikes a 
hard quark which carries a momentum fraction of the Lorentz contracted proton with large positive 
light-cone momentum P+. The quark is turned around and then exits the proton. Being considerably 
off-shell from the hard collision, it begins to shower gluons. The initial radiations have large transverse 
momentum l± (due to the large virtuality), and thus have shorter formation times: 



Where y is the momentum fraction of the parent parton carried away by the radiated gluon. Later 
radiations (those with larger formation times) have smaller transverse momentum. While not illustrated 
in the figure, the radiated gluons are also off-shell and tend to radiate gluons with smaller off-shellness. 
Eventually the local virtualities of the partons is so low that one cannot apply a perturbative partonic 
picture and hadronization begins to set in. As a result, the collection of hard collinear partons turns into 
a coUimated jet of hadrons. The leading hadron in this picture is the one with the largest longitudinal 
momentum and, given some form of localized momentum conservation, is the result of hadronization of 
the highest longitudinal momentum part or the largest formation time part of the jet. 

The picture in p-p is somewhat similar except that one has at least 2 back-to-back jets being 
produced. Both of these jets are produced considerably off-shell and tend to lose this virtuality by 
successive emissions. Unlike the case of DIS, jets which arise from the produced hard partons in large 

p-p collisions have a large momentum transverse to the incoming hard partons. We will always denote 
the momentum in this direction as Pt, differentiating it from the momentum of gluons transverse to 
the produced jet which will be indicated with a _L, as in l±. The factorized formula to obtain the cross 
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Figure 1: DIS on a nucleon leading to the formation of a quark jet which showers gluons in 
vacuum leading to the formation of a coUimated jet of hadrons. 



section of single hadrons is given by the first term on the right hand side of Eq. 
this formula simplifies to, 



At leading order, 



dyd?pi 



dXadXb GaiXa, Q)Gb(Xb, Q) 



ddab^cxiQ , XaP, XbP, Pt/ Z) D^{z,Q) 



dt 



TXZ 



(9) 



The reader will note that we have chosen all factorization and renormalization scales to be Q which 
is the hard scale in the problem. At leading order all hard scales that appear in the calculation are 
equivalent; in most cases one simply picks Q — Pt- 

We now consider the process of DIS on a large nucleus. This will often be referred to as A-DIS (A 
referring to a nucleus). In the Breit frame, one may neglect the soft interactions between the various 
nucleons as these occur over a long time scale. The nucleus may then be modeled as a weakly interacting 
gas of nucleons traveling in the positive light-cone direction. The virtual photon strikes a hard quark 
in one of the nucleons and turns it backwards as in the case of the DIS on a single proton. The quark 
then propagates through the nucleons directly behind the one that is struck, as illustrated in Fig. |2l 
where we draw the struck quark as propagating outside the line of nucleons for clarity. The quark is 
virtual at production and radiates a shower of gluons with progressively longer formation times, similar 
to the case of DIS on a proton. In this case however, both the quark and the radiated gluons tend to 
scatter off the soft gluon field in the nucleons. This is indicated by the zig-zag lines. Note, the zig-zag 
lines do not indicate a pomeron or double gluon exchange, but rather are single-gluon exchanges which 
are distinct from the gluons in the shower of the jet. Diffractive exchanges with the nucleons may be 
neglected in the case of very large nuclei. This will be justified in the next chapter. 

The multiple scattering of the partons in the shower changes their momentum distributions and 
as a result the final hadronization pattern is modified. In the case of a jet produced in a heavy 
ion collision, the picture is qualitatively similar except for the absence of nucleons. The different 
components of the jet now scatter off the quark-gluon substructure of the degrees of freedom in the hot 
deconfined matter. The leading effect of the multiple scattering on the shower profile is a broadening 
of the distribution in transverse momentum (recently this has been experimentally observed in cold 
nuclear matter [23]). Besides transferring transverse momentum, the medium also exchanges energy 
and longitudinal momentum with the jet (in the language of light-cone momentum, both (+) and (— ) 
components are exchanged between the jet and the medium). While energy and momentum exchange 
(including transverse momentum) in the right proportion may cause minimal change in the virtuality 
of a given jet parton, arbitrary momentum exchanges may noticeably change the virtuality leading to 
induced radiation. In the case that the energy (and longitudinal momentum) of the jet partons far 
exceeds that of any constituent in the medium, such exchanges lead to a depletion of the light-cone 
momentum in the forward part of the jet. This is often referred to as radiative energy loss and results in a 
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Figure 2: DIS on a large nucleus. A quark in a nucleon on the front side is struck and 
then propagates through the nucleons behind the struck nucleon. In the process both the 
quark and the ensuing shower gluons scatter off the soft glue field inside the nucleons. This 
modifies the shower pattern. 



suppressed yield of leading particles [211 ISS |26l [271 1211 |29l [30] . This suppression has been quantitatively 
established in both DIS on a large nucleus [31] and in jets in heavy ion collisions [21 [3]. The primary 
difference between the two cases is the fact that a heavy-ion collision is not a static environment and 
evolves rapidly with time and then finally disintegrates into a cascade of hadrons. 

The calculation of the modification of the hadronization pattern due to the multiple scattering in- 
duced radiation will be the subject of the next two chapters. Until recently, experiments have only 
been able to measure the modification of the yield of the leading hadron or the correlation between 
the leading and next-to-leading hadrons. The basic theoretical object required to describe this aspect 
of jet modification is called the medium modified fragmentation function. Along with the perturbative 
calculation of the evolution of the fragmentation function in vacuum, one computes the change in the 
fragmentation function due to the broadening and stimulated emission that occurs as the jet passes 
through matter. The medium modified fragmentation function is now a function of not just the mo- 
mentum fraction and the scale Q^, but also of the energy of the jet and the distance travelled in the 
medium. While it has become conventional to call this a "medium modified" fragmentation function, 
it should be pointed out that in all pQCD based calculations, the fragmentation is assumed to occur 
once the jet has escaped the medium. 

We should point out that while this review will describe jet modification in dense matter as an 
extension of the factorized processes in vacuum pQCD [this is often referred to as the Higher-Twist (HT) 
scheme [32 1 [33 1 [M]]. there are other approaches to this problem. An entirely orthogonal approach is that 
based on finite temperature field theory based on the work of Arnold, Moore, Yaffe and collaborators [33 
[SSI [33 [3S1 [33 [ID]- Referred to as the AMY scheme, this formalism considers the effective Hard Thermal 
Loop (HTL) effective theory of dense matter and considers the hard jet to have the same virtuality 
or mass scale as a thermal plasma particle but with energy E ^ T. Following this, one identifies 
and resums the collinear enhanced contributions emanating from the scattering and induced radiation 
off the hard parton. The entire calculation is carried out at the scale of the temperature T which is 
assumed to be large (T — )■ oo) so that the effective coupling g(T) is small [giT) — >■ 0]. Another approach 
developed by Armesto, Salgado, Wiedemann and collaborators models the medium as a series of Debye- 
screened, heavy, colored scattering centers. In this approach, referred to as the ASW scheme, the hard 
parton radiates a virtual gluon which is then progressively brought on shell by a large number of soft 
scatterings off these heavy centers [HI [HI [131 [HI [H] . Yet another approach developed by Gyulassy, 
Levai, Vitev and their collaborators considers the same medium as ASW, however both the hard parton 
and the radiated gluon undergo a few but hard interactions with the centers leading to the emission of 
the gluon [l6l W7\ [HI [HI [501 [SI] • After the calculation of the single gluon emission kernel, the AMY 
scheme uses rate equations to incorporate multiple emissions whereas the GLV and the ASW use a 
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Poisson emission Ansatz. 

The comparison of the yield of leading particles between the case with and without a medium is a 
measure of the properties of the medium as felt by the jet. In DIS, one measures the fragmentation 
function to produce a hadron with a momentum fraction z of the original quark momentum immediately 
after being struck. For the case where the virtual photon carries a momentum, q = [— (5^/(2g~), g~, 0] 
and the incoming struck quark has p = [Q"^ /2q~ ,0,0], the outgoing quark has a momentum [0,g~,0]. 
Thus the hadron momentum is ph — zq~ . The fragmentation function for DIS on a proton (nucleus A) 
and the ratio R are defined as, 

( \ ^ ^^P(^) A G ap daA/dz 

Dp(^A){z) = 1 — , and R= — - — — . (10) 

(Jpi^A) dz a A dap/dz 

In the equation above, ap^A) is the differential DIS cross section off a proton p (nucleus A) integrated 
over a limited range of and energy imparted from the electron v = q~ . In order to compare with 
theoretical calculations, the fragmentation function in the case of DIS on a large nucleus Da{z) should 
be identified with the medium modified fragmentation function D{z). 

In the case of p-p or heavy-ion collisions, the jets span a range of momentum depending on the 
momenta of the two partons that undergo the hard scattering. In this case one cannot isolate separate 
bins in the final state momentum fraction z and thus cannot measure the fragmentation function directly. 
Instead one measures the ratio of the binary scaled differential yield to produce a high p^ hadron in a 
heavy- ion collision to that in a p-p collision referred to as the nuclear modification factor Raa- The Raa 
may be measured both differentially as a function of angle with the reaction plane and the transverse 
distance between the centers of the colliding nuclei, the impact parameter b, or integrated. The angle 
integrated Raa [also integrated over a small range of impact parameter (femm to bmax)], is defined as 



min I'-'max } 



where, (A^bin(^)) is the mean number of binary nucleon nucleon encounters per ion-ion collision in the 
range of impact parameters chosen. The invariant yield in a heavy-ion collision and in a p-p collision 
are related to the the invariant differential cross section by the relation, dN = da /a, where the total 
cross section for a heavy-ion collision may be estimated as the geometrical cross section in that range of 
impact parameter selected. The relation between A'^bin and the nuclear density profile will be discussed 
in Chapter [5l 

The calculation of the nuclear cross section to produce a hard hadron may be expressed as the 
convolution of the nuclear PDFs, the hard partonic cross sections and the medium modified fragmen- 
tation function. The jets may be produced at various locations in the hot matter and propagate in any 
direction in the transverse plane. The location of the production point and the direction determine 
the extent and intensity of the medium as felt by the jet. One thus needs to integrate over all allowed 
production points and directions. This procedure will be described in Chapter |3 



1.4 Medium transport properties: what can be learnt from jets 



As described in the preceding subsections, a jet is essentially a collection of high momentum particles 
which are somewhat coUimated in a given direction. In the partonic part of the jet, these are virtual 
partons. To make this explicit we ascribe a virtuality /i^ ^ to the hard jet, where describes the 
scale of the hard interaction which produces the jet. In effective field theory approaches to hard jets 
such as that of Soft Collinear Effective Theory (SCET) [521 [53] one introduces a small parameter A with 
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Proton 



Figure 3: A single Coulomb (or Glauber) gluon interaction between a hard jet and a proton 
in the nucleus. 

A ^ 1; colhnear radiations have a transverse momentum of ~ \Q and the virtuality of the jet may 
then be surmised as \^Q^. Thus /i^ ~ {\QY\ the two terminologies will often be used interchangeably 
in this review. 

As the jet passes through matter each of these particles will scatter off the various constituents 
that it encounters. In the Breit frame (or infinite momentum frame) the nucleus has a large boost in 
the direction opposed to that of the jet. As we pointed out earlier, one may consider each nucleon 
within the nucleus to be traveling in almost straight lines independent of each other. At sufficiently 
high energies, the "large-x" partons in the nucleons may also be considered to be traveling in straight 
lines independently of each other (see Fig. |3]). The hard virtual partons in the jet with virtuality /i^ 
will resolve this sub-structure down to transverse sizes of order 1///^ ~ (AQ)"^. The hard partons in 
the jet will scatter of these partons by exchanging gluons with transverse momenta /cj, ~ ~ \Q. As 
a result, the hard partons in the jet will undergo a transverse diffusion as they propagate through the 
extended matter. 

Along with the exchanged transverse momentum there may also be a certain amount of negative 
light-cone momentum {k~) which may be exchanged between the partons in the jet and those from 
the medium. In this review, we will refer to this as "elastic energy loss" or "drag", even though the 
incoming parton from the medium is not on-shell, is not a quasi-particle of the medium and may not 
go back on shell after the scattering. This is primarily done to distinguish this type of light-cone 
momentum loss from the light-cone momentum loss that occurs due to radiation. This is somewhat 
different from what is traditionally referred to as elastic energy loss, which involves a certain energy 
and momentum transfer to an in-medium quasi-particle which remains as one quasi-particle after the 
interaction, but with larger energy. The remaining component of the transferred momentum is 
completely determined by insisting on the criterion that the jet parton virtuality does not change too 
much after the momentum transfer. Thus if the original off-shellness is 2q^q~ ~ /i^, then after the 
transfer, one obtains the virtuality as, 

2(g+ + fc+)(g- + k-) - fc2 = + 2^k- + fc+(g- + fc") - kl. (12) 

Q 

Thus, for k~ <^ we obtain that the virtuality of the parton will remain more or less unchanged if 
/c"*" = k^/^l^q" + k^)]. For cases where /c"*" exceeds this value, the hard parton from the jet is taken 
further off-shell leading to an induced radiation. 

Treating the case of momentum change of the hard parton due to induced radiation separately, we 
obtain three separate components of momentum being exchanged between the jet and the medium: two 
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components of transverse momentum and one component of negative light-cone momentum. At every 
exchange there is a distribution in transverse and hght-cone momentum being imparted between the 
jet and the medium. Given the large number of exchanges, we invoke the Gaussian approximation, i.e., 
we approximate these distributions to be Gaussian and consider only the mean and the variance. The 
Gaussian approximation based on the central limit theorem is not completely unjustified. Each parton 
interacts multiple times with the medium and the shower distributions contain several hard partons. 
Along with this, one should also consider that except for the highest pt (or highest energy) hadrons, 
each bin in hadron momentum (or momentum fraction) contains several events that need to be summed 
over. 

Given cylindrical symmetry around the jet axis, on may further argue that the imparted kj]_ = ky'^^, 
and the mean of the transverse momentum distribution is vanishing. Thus the first transport coefficient 
may be defined as the variance of the distribution of imparted transverse momentum per unit length 
traversed. This is referred to as q, 



l^x I I r + 1^' 



y±\L 



L 



(13) 



Where is the total transverse momentum gained in traversing a length L. In the case of negative 

light-cone momentum exchange, the mean of the distribution yields the drag per unit length (e), whereas 
the variance is the fiuctuation in light cone momentum transfer per unit length (62): 

e = — , and 62 = ^ (14) 

Given these two quantities, the Gaussian distribution of the negative light-cone momentum transfer is 
completely specified. 

Thus the modification of hard jets will reveal at most these three quantities. Due to the large number 
of scatterings that need to be summed over in order to obtain measurements with small error bars, 
access to higher moments in these distributions is somewhat limited. These transport coefficients can 
in principle be calculated given a model of the medium, e.g., in the asymptotically high temperature 
limit, the system can be described in the weakly coupled quasi-particle picture afforded by leading 
order Hard Thermal Loop (HTL) effective theory. This derivation will be described in some detail in 
the remaining review. Even in this case, the value of the transport coefficients are unknown unless the 
value of the coupling is specified. The relevant coupling is the in-medium coupling between thermal 
quasi-particles which is a priori unknown and is therefore set by fitting calculations to at least one data 
point. It should be pointed out that if the medium is determined to not be weakly coupled, the HTL 
formulae cannot be applied and some other model will have to be used. The value of the transport 
coefficient, in principle, will have to be evolved up to the scale of the hard jet prior to use in a jet 
quenching calculation. 



1.5 Other approaches and medium response observables 

In this review, we will primarily address the question of how a hard jet can be used as a weakly 
coupled probe to study the properties of a dense QCD medium. As such we will solely be interested 
in how different properties of the medium, codified as a set of transport coefficients will modify the 
shower pattern of the jet. As mentioned above, throughout this review, we will solely describe the 
the modification of the partonic portion of the jet shower, i.e., that part which may be computed 
using pQCD. This is often referred to in the literature as partonic energy loss to distinguish it from 
an alternate mechanism of hadronic energy loss. In this alternate hadronic energy loss scenario, one 
assumes that the jet hadronizes within a very short distance (of the order of 1 fm) and produces a 
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shower of hadrons. These hadrons then muhiply scatter and lose forward momentum in the dense 
medium which is also assumed to be hadronic [SU 

When applied to heavy-ion collisions, these calculations assume a minimal partonic phase, if any. 
While such theories have experienced a degree of success in the description of jets produced in DIS on 
a large nucleus [31], they have performed less favorably in comparison to RHIC data [56]. A slightly 
more successful variant has been the model based on pre-hadronic absorption [571 EH [59] which tend 
to model the absorption of the QCD string prior to hadronization. No doubt, all three mechanisms 
will be present in any description of jet modification where one attempts to capture the fate of all the 
hadrons which materialize from the shower, in media of arbitrary length. However, in sufficiently short 
media, and for considerably energetic jets, one may assume that the leading hadrons have hadronized 
outside the medium (these statements will be made more quantitative in the ensuing chapters). The 
yield of these leading hadrons may then be described by a factorized vacuum fragmentation function 
convoluted with the medium modified distribution of hard partons which have escaped the medium. 
The modification of the distribution of these few hard partons may then be calculated using pQCD 
with a minimal amount of non-perturbative input in the form of transport coefficients. The vacuum 
fragmentation functions are obtained from e^e~experiments and are thus well known. Thus the study of 
this subset of observables involves only one set of unknowns: the transport coefficients of the medium. 
It is in this very restricted sense that pQCD based jet modification may be used as a probe of the 
medium. This review will be focused on this very restricted set of observables. 

As the medium modifies the jet, the jet in turn modifies the medium. There is a growing consensus 
that part of the energy radiated by a jet is deposited in the medium and this modifies the evolution of 
the soft medium at that location [601 [61]. Given the "supersonic" velocity of the hard jet, this may lead 
to the production of a Mach Cone [52] which trails the hard jet and may be observable as an excess in 
the hadronic yield at a large angle to the associated away side jet in the final soft hadronic spectra that 
is associated with a hard jet trigger [631 [M]- 

While the response of a soft medium to a hard jet may not be completely calculable in pQCD, in 
recent years an alternate theory based on the Anti-DeSitter space Conformal Field Theory (AdS/CFT) 
conjecture [55] has been used to compute both the drag experienced by a heavy quark [BSl [SZ] and the 
Mach cone left in a wake of of a hard jet [HH]- Unlike calculations based on pQCD, these theories 
assume that the hard jet is strongly coupled with the dense medium. In the absence of a pQCD 
based picture it is not clear if the final hadronization process or the initial production process may be 
consistently factorized from the energy loss calculation. The three topics of (pre-) hadronic energy 
loss, medium response to a hard jet and the alternate theories of energy loss based on the AdS/CFT 
conjecture will not be covered in this review. For the reader interested in such topics we recommend 
these other excellent reviews: Refs. [701 [HI [721 [73] . 
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2 Scattering induced single gluon radiation. 



In order to compute the energy loss of hard jets and the medium modified fragmentation function, one 
needs to compute a series of multiple scattering and multiple emission diagrams. In all formalisms, the 
means to achieve this is to first compute the single gluon emission kernel due to multiple scattering. This 
is then iterated to include the effect of multiple emissions. The methodology of iteration is somewhat 
related to the approximations made in the underlying calculation of the single gluon emission kernel. 
In this chapter, we will describe this calculation in some detail. The iteration of the kernel will be dealt 
with in the next chapter. As stated in the introduction, our description of the single gluon emission 
kernel will follow that of the higher-twist approach. The differences with the other approaches will be 
pointed out at the end of this chapter. 

We present the formalism in the Breit frame. In the case of DIS on a large nucleus this is char- 
acterized by the frame where the absolute magnitude of the (-I-) and (— ) components of the photon 
momentum q and the large (-I-) component of the initial momentum p of the struck quark are equal: 

g = 1^--^, -^,0^ and p+ = a;BP+ = ~ ^ 0. (15) 

is the momentum of a proton. The struck quark has a final (— ) momentum ~ = Q/\/2 and 
travels in the negative z direction. It now scatters off the gluons in the nucleons which follow the struck 
nucleon. The picture is similar to that in Fig. |31 

We will present the calculations in negative light-cone gauge A~ = 0. This makes the discussion 
of gluons emitted coUinear to the outgoing quark particularly simple. The first step is to quantify the 
magnitude of the different components of the gluons being exchanged between the outgoing quark and 
the incoming nucleons. The final outgoing quark has a negative light-cone momentum of qj ~ Q. Being 
close to on-shell, it has a virtuality qj ~ A^Q^, which is built up from some combination of a gj' ~ X^Q 
and a small transverse momentum qf± ~ \Q. 



2.1 DIS on a proton and vacuum radiation 



To familiarize the reader with the notation, we compute the simple process of vacuum gluon radiation 
from a quark struck by a hard photon. This constitutes the primary contribution to the scale evolution of 
fragmentation functions in vacuum. With this process we will also demonstrate an elementary example 
of factorization at leading order, where we factorize the parton distribution function (PDF) from the 
hard cross section. The basic quantity to be computed is the differential cross section for an electron with 
an initial momentum Li to scatter off a proton (with momentum P~^) with final electron momentum 
L2 and producing an outgoing quark with momentum Ig = [|Zgj_p/{2g~(l — y)},q~{l — y),lqj_] and a 
gluon with momentum / = [\l ±1"^ / {2q~ y) , q~ y , I ±]: 

£(Li) + p(P+) C{L2) + Q(/,J + G{h) + X. (16) 

The complete cross section for the process is represented by Fig. HI The rectangular blob at the bottom 
of the diagram represents the contents of the proton after the hard quark q is struck by the photon and 
removed from the proton. 

The matrix element that needs to be evaluated may be represented as 

M = {X,l,lg;L2\Texp 



dtHi{t) 



\P;Li). 



;i7) 
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Figure 4: DIS on a nucleon leading to the formation of a quark and a radiated gluon. Left 
panel represents the complete diagram whereas the right panel represents the hadronic tensor. 
The dangling photon lines do not represent photon propagators, but rather the scattering of 
the quark off a photon. The photon propagator is not contained in the hadronic tensor, see 
text for details. 



The state |X) represents an arbitrary hadronic state. For the case of interest, \L2,l,lq, X) represents a 
state with at least a quark and a gluon endowed with the requisite momenta and an outgoing electron 
with momentum L2. At lowest order \X) is simply the vacuum state. In the equation above, Hj{t) 
represents the interaction Hamiltonian in the interaction picture for QCD and QED. The eventual 
evaluation of the exponent will involve two orders of the electro-magnetic (EM) interaction and all 
orders of the strong interaction. For the case of one gluon emission in this section we will only have to 
expand the exponent to one order of the strong interaction (and two orders of the EM interaction) i.e.. 



00 POO 



M = (X,/,/g;L2|^ / d%llJeiye)eYAf^{ye)lpeiye) I c/^yo^(yo)e7''^M(?/o)^(?/o) 



00 



'3! 

X / d^y^^P{y^)gYt''At{y,y4^{y,)\P;L,). 



In the equation above, '0e is the wave function of the lepton, while '0 is that for a quark. The restriction 
to only two orders of the EM interaction, constrains the calculation to the one photon approximation. 
Squaring the matrix element one may calculate the cross section for this process. It may be straight- 
forwardly demonstrated, in the one photon limit, that the process can be decomposed into a purely 
leptonic part and a partonic (or hadronic) part connected by a single photon propagator. The leptonic 
part will involve the trace over the Dirac matrix structure of the electron momenta; we will simply write 
this down without proof. The hadronic part will be somewhat more complicated as it will contain both 
a perturbative component describing the scattering of a quark off the virtual photon followed by the 
emission of a hard gluon and a non-perturbative component which describes the probability to find a 
quark with a particular momentum fraction in the incoming nucleon. This component will be described 
in some detail. 

The quadruple differential cross section for the process in Eq. f|T6|) . in terms of L2,l±,lqj_,y can 
be decomposed, as stated above, into a purely leptonic part and a hadronic part in the one photon 
exchange approximation. This is symbolically represented as. 



d^L2dHqJ^±dy ^'"^^^'^ dHq^dH^dy ^^^^ 

where L"^ and W^^ are the leptonic and hadronic tensors respectively. The Mandelstam variable 
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s = {p + LiY. Each of the factors of Qaii/Q'^ and g^vlQ'^ represents a photon propagator, one from the 
amphtude and one from the complex conjugate. 

The leptonic tensor has the straightforward definition, L"^ = Tr[/^i7" /^27^]/2, while the hadronic 
tensor, at leading order in the strong coupling, may be decomposed as 

Wr = j Ao(P|^(yo)7^^7XO)|P) = j ci^i/oTr[l-7'^^7li^(yo)0°°(yo). (20) 

This expression requires some explanation. The interaction terms in Hj are those contained in the QCD 
[Eq. ([1])] and QED Lagrangians. These are composed of solely quarks, gluons, leptons and photons. 
While we have, somewhat artificially, constrained the final state to be a quark and a gluon, the initial 
state is a proton. In the Breit frame, the high energy proton can be approximated as a weakly interacting 
gas of coUinear quarks and gluons. The distribution of these partons in bins of light-cone momentum 
fraction x depends on the transverse size of these partons (i.e., resolution of the probe). In this high 
energy limit, the projection of the product of ip\X) and its complex conjugate {{X\%1)) along the large 
momentum direction may be viewed as an annihilation and creation of a near on-shell quark in the 
proton's wave function. The Fourier transformation of the expectation of this operator product is 
referred to as the parton distribution function (PDF), 

f{x) = J dy^e~^^^^yoF{y,) = J rf|/-e-^^^o-(p|^(^-, 0)1^^(0, 0)|P). (21) 



Unlike the parton distribution function in Eq. (12T1) above, the function F{y^) that appears in Eq. (1201) 
is not yet on the light-cone. Note that in the equation above only appears. This happens after the 
invocation of the high energy or collinear approximation. The incoming parton is assumed to be endowed 
with very high forward momentum {po^ = xqP~^,Pq — >■ 0) with negligible transverse momentum po± 
Pq. Within the kinematics chosen, the incoming virtual photon also has no transverse momentum. As 
a result, the produced final state parton also has a vanishingly small transverse momentum {i.e., with 
a distribution that may be approximated as S'^{p±) up to corrections suppressed by powers of the hard 
scale). In this limit, the leading spin projection of the pieces which represent the initial state and final 
state may be taken. The factors, 

are used to obtain the spin projections along the leading momenta of the outgoing state and the incoming 
state. See Refs. [SU [71] for precise details of how this is done. 

The other function in Eq. fl20|) is O^^ which represents the physics of the final state after the hard 
scattering with the photon. The superscript on the operator O^^ implies that the quark undergoes no 
scattering in the final state. Taking the leading projection in Dirac matrix structure we obtain 



000 ^ rj^ 



2_0OO 



(27r)4 (27r)4 (27r)4 (27r)4 



2 [po + q)-^ - le ^ {pq + qy + le 



Tr[ef]. (23) 



In the equation above, Ga/sil) is the sum over the product of polarization vectors of the radiated gluon 
with momentum I. In light cone gauge A ■ n = A~ = (with n = (1, 0, 0, 0)) this is given as, 

Gai3{l) = -g^p + -j-^ • (24) 
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To evaluate the above Feynman integral one integrates over the internal locations z', which will 
yield the momentum conserving delta functions. This will set po = p'o- Then one approximates the 
numerators of the fermion propagators as 'J^q~ ■ The denominators may be expressed as 2P~^q~{xQ — 
XB^i^) where xq = p'^/P^. The sets = l\/2l~ . The other 5-function over Iq^ sets xq = xb + xl 
where xl = l'^/[2P^q~y{l — y)]. The remaining steps involve carrying out the simplifications associated 
with the 7 matrices and carrying out the integrations over the momenta that do not appear in the 
integrand. This yields 

O'' = 5(y^)5^(yoJe-(— )-^^o-^ / dydll2-2y + y^ _ ^^^^ 

271 J lj_ y 

The resulting phase factors constrain F{y^) to the light cone and convert it into the PDF. The factor 
{2 — 2y + y'^)/y is the splitting function P{y) which expresses the probability for the hard quark to 
radiate a gluon. Re-incorporating this expression back into Eq. (120 p . we obtain the differential hadronic 
tensor for a quark in a proton, struck by a hard virtual photon, to radiate a gluon with momentum 
[q~y, l±] and itself have a momentum [q~{l — y), lq_^]: 



iyi,lin,, - T. ^<M^^ + ^.)i^Piy)^%. + rj. (26) 

The above result should be compared with the simpler result of DIS without any radiation in the final 
state. The derivation is similar to that presented for the case of a produced radiation, we here simply 
state the result (referring the interested reader to Ref. [33] )• The differential hadronic tensor for a quark 
in a proton, struck by a hard photon, to have a momentum [q~, Ig^] is 

^ = J2^7rQlfq{xB)S%^). (27) 

a Iq^ ^ 

In the equation above, fq{x) is the quark PDF defined in Eq. (1211) . Note that the arguments x in 
Eqs. (l26|27p are slightly different. They differ by the quantity xl = l']_/\2P^q~y{l — y)]. This should be 
understood as the ratio of the off-shellness of the outgoing quark to Q"^ /xb = s, i.e., 

XLS = — = XL2P^q- = . = 2(/+ + /+)(/- + /-) - /"I + lq^\ (28) 

Xb y{l-y) " 

The momentum component p'^ = x^P^ is the extra (+)-component of the momentum that must be 
brought in by the incoming quark such that the quark after being struck should be off-shell enough to 
radiate the gluon with momentum \L\/{2q~y),q~y,l±\. We work in the limit where a collinear jet is 
produced in the hard interaction, so that the off-shellness is very small compared to the hard scale Q^: 

XL2P+q~ = < = XB2P+q- =^ xl < xb- (29) 

As mentioned in the introduction, this ratio of the off-shellness of the produced quark to is used to 
qualify the small parameter A, i.e., ~ AQ. In this limit of transverse momenta, xl = X'^xb{^xb) if 
y ~ 1(A). In either case, this small correction to in Eq. fl26|) can be neglected. 

In both equations fl26|27p . the distribution of the final outgoing quark involves a 5-function. This is 
meant to be a simplification. A better expression is to use a narrow, normalized Gaussian distribution. 
The width is of order ~ XQ, i.e., small compared to the hard scale but still perturbatively large. The 
width of this outgoing quark distribution is related to the transverse momentum and off-shellness of 
the incoming quark; this controls the virtuality of outgoing quark and thus the virtualities involved in 
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all final state processes. This also controls the scale at which the coupling constant is evaluated in the 
case of final state gluon emissions. 

In the limit of small xl-, the ratio of the cross sections for both processes, computed using Eq. f|T9|) . 
assumes a rather simple form. Integrating out the transverse momentum of the produced quark gives 
the differential number of gluons radiated, 

dNg ^ 1 c?er[g+^^g+g+x] ^ (^s{1]_)Cf p^y^ 
dydl\ cr[q+^^g+x] dydt]_ 2-kI']_ 

In the subsequent section, we will derive the cross section for a hard quark to radiate a single gluon 
while undergoing multiple scattering in the medium. This will involve modifying a related version of 

We close with a comment regarding Eq. (130|) . The reader will note that the expression for the 
differential number of gluons in y, l± possesses both an infra-red y ^ [P{y) 2/y as y ^ 0] and a 
coUinear singularity l± — )■ 0. In the calculation of the final hadronic distribution, one convolutes this 
gluon distribution and the related final quark distribution with a fragmentation function D as shown 
below in Eq. ( 1311) . In this convolution one also has to include the effect of virtual corrections, i.e., 
instances where the leading parton radiated a gluon and then reabsorbed it. Such diagrams also contain 
an infrared and collinear divergence. In the convoluted expression to obtain the yield of hadrons, the 
infrared singularity cancels between real and virtual diagrams. The real splitting function is replaced 
with the well known (+)-functions [18j. The collinear divergence, however remains even after the 
inclusion of virtual corrections. The source of this divergence are gluons with lj_ — )■ 0, and thus those 
with formation times Tf — oo. In a factorized approach, such long distance effects are absorbed into 
a renormalization of the final fragmentation function to produce hadrons. In the case of single hadron 
inclusive DIS, absorbing gluons with transverse momenta up to the final factorization scale l± = fi <^ Q 
yields the scale dependent fragmentation function [D{z,n'^)] to produce a hadron with momentum 
fraction z = Ph/Q- Including the effect of multiple emissions yields the DGLAP evolution equation for 
the fragmentation function: 



91og(/i2) " 27r 



Jz y \y J 



(31) 



In the in- medium version of Eq. fl5U]) . Landau-Pomeranchuck-Migdal (LPM) interference tends to 
cancel the collinear divergence making the differential yield of gluons finite at — 0. As a result, 
one obtains a finite energy of gluons emitted on integrating over Two of the formalisms use this 
value to construct an iterative formalism, thus ignoring the fact that a fraction of the number of gluons 
are produced much later in time and should not be included in the calculation. This will be further 
discussed in Chapter [3l 



2.2 Multiple scattering induced single gluon radiation 



In this section, we compute the cross section for single gluon emission induced by multiple scattering. 
Iterating this process and convoluting with the fragmentation function will produce the medium mod- 
ified fragmentation function. We will approximate each scattering to only transfer a small transverse 
momentum k± ~ XQ and a (+)-component ~ X^Q between the jet and the medium partons. There can 
also be an exchange of a small (— )-component; this leads to elastic energy loss and will be discussed 
separately in Chapter HI 
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2.2.1 Multiple scattering without gluon radiation 



We begin by considering the case of a hard quark produced in a DIS on a large nucleus, which 
produces a hard quark propagating through the nucleus encountering multiple scattering without ra- 
diation. The virtual photon strikes a hard quark in one of the nucleons and turns it back towards 
the nucleus. The quark then propagates through the nucleus without radiating. In this process the 
hard quark scatters n times in the amplitude and in the complex conjugate. This is represented by 
the diagram in the left panel of Fig. |S1 The diagram represents a quark in one of the incoming nu- 
cleons with momentum pg, being struck with the photon with momentum q = [—XbP~^, 0, 0]. The 
outgoing quark has momentum q[ = p'q + q. After encountering the j*^ scattering, its momentum is 
q'j+i =Po + q + Ei=i K- 

In the case of the higher-twist scheme, these gluon lines represent the gluon field at the point at 
which scattering takes place. Thus there is no meaning to crossing of gluon lines. The sources of the 
gluon lines are the nucleons in the nucleus (or rather the partons within these nucleons). The only 
assumption made regarding the gluon momenta are simple scaling relations regarding the momentum 
of the exchanged gluons and the magnitude of the corresponding components of the vector potentials. 
This is somewhat different in the other schemes. Under the assumption that the gluons transfer a 
transverse momentum of k±_ ~ XQ, and the incoming and outgoing quark lines remain close to on-shell, 
the component is constrained with the leading contribution given by the equation, 

(q + ky = 2q~k+ ~kl = =^ k+ = -^r^\'^Q. (32) 

2q- 

Ascribing the same criteria to the hard quark or gluon in the nucleons from which the exchanged gluon 
originates, implies that k~ ~ k\/2p'^ (where p'^ here is the momentum of some hard parton) also scales 
as ~ X^Q. Exchanged gluons with transverse momenta much larger than their longitudinal components 
are referred to as Glauber gluons or Coulomb gluons. If the (-l-)-component were to become larger than 
kj_/2q~, it would drive the jet parton to go off-shell and radiate a hard gluon. If the k~ component 
were to become larger, this will lead to energy loss from the jet parton but the parton in the nucleon 
will then go off shell and may radiate. As jet quenching measurements do not concern themselves with 
the fate of the target, the latter sort of momentum transfer is often referred to as elastic energy loss for 
the jet. 

Given these approximations we may now compute the scaling of the different components of the 
gauge field. To remind the reader we are calculating in the Breit frame where the hard quark moves in 
the —z direction with a large light cone momentum q~ ^ Q and the valence quarks inside the nucleon 
are moving in the +z direction with a large ~ Q. We use the linear response formula to ascertain 
the power counting of the A"*" field. Suppressing the color superscripts we obtain, 

A^{x,) = j d'y,V'^''ix,-y,).My,). (33) 

In the equation above, V is the gluon propagator and at leading order in the light cone gauge A" = 
is given as. 

In Eq. f l33|) . J^{yi) = 4'{yi)l'^'^{yi) is the current of partons in the target which generates the gluon 
field. The fermionic operator may be decomposed as, 

Hvi) = I ^'iPKe-^'-'' + ■■■ (35) 
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Figure 5: Left panel: DIS on a nucleus leading to the formation of a quark jet which which 
is constrained to propagate through the nucleus without radiating. Its only interactions are 
scatterings. Right panel: The dominant length enhancement arises from nested scattering 
diagrams, where the blobs represent individual nucleons, see text for details. 



The scaling of the fermionic operator depends on the range of momentum which are selected from the in- 
state by the annihilation operator. Note that this influences both the scaling of the annihilation operator 
Qp as well as the bispinor u{p). The power counting of the annihilation operator may be surmised from 
the standard anti-commutation relation, {ap,^^/^} = (27r)^(5^(p — The power counting of the 

bispinor can be obtained from the completeness relation, 'Yls'^p^p ='t> = 1~V^ + 1^V~ ~ 7± ' P±- 
Substituting the equation for the current in Eq. ( 133|) . and integrating out y, we obtain. 

If the incoming and out going momenta p and q scale as collinear momenta in the (+)-direction, i.e., 
p ~ Q(1,A^,A), then we get, — j?) ~ [A^Q'^]"^, as one of the momenta will involve the large 
(+)-component and the remaining are the small transverse components. Thus the annihilation (and 
creation) operator scales as \~^Q~^^'^ . Also in the spin sum ~ Q ^^id thus u{j)) ~ ~ Q^^^; one 
can check that the 7"*" projects out the large components in u and u in the expression u{q)'y~^u{p). We 
also institute the Glauber condition that p+ — g+ ~ A^Q, p~ — q~ and p± — q± ^ XQ. 

Using these scaling relations we correctly find that the bispinor scales as XQ^^^. However, to obtain 
the correct scaling of the gauge field one needs to institute the approximation that q^ = p'^ + 
where k'^ ~ X^Q. This condition is introduced by insisting that the (+) momentum of the incoming 
and outgoing state, which control the scaling of and a^, are separated by k'^ ~ X^Q. This is used to 
shift the dq^ — >■ dk^ and as a result we obtain the scaling of the field as X^Q. Following a similar 
derivation, with the replacement 7+ — )■ 7_l we obtain the scaling of the A± ~ X^Q. Thus in the Breit 
frame and A~ = gauge, the A'^ field is dominant over the A± components which may be neglected. 

The evaluation of the right panel of Fig. |5] is straightforward within the approximations outlined 
above. The numerator of every propagator is replaced with 7^g~ and every vertex with 7"/!+. The 
denominators are simplified by contour integration, e.g., for the simple case of the first propagator, we 
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would get, 



271 {q- + hy " 



27r q-k+ - {kiy + ie 



(37) 



Thus carrying out the contour integrations, simply orders the locations of the scatterings and sets the 
(+) components of the momentum to be k']_/2q~ . Carrying out these integrations on both sides of the 
cut we obtain an expression which is a modification of Eq. fl?l 
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In the equation above, we use the short hand 



E 2pi ■ 

,fc=0 



/2P+g to save writing, and 



similarly for the complex conjugate. 

Further evaluation requires two sets of approximations. The first arises from the large length limit 
of a large nucleus. Each of the dy~ or dy'f integrals may be extended up to the light-cone length of 
the nucleus L~ . The largest factor of L~ is obtained in the limit that the gluon operators A~^{yi) are 
distributed over the largest number of nucleons. The color singlet nature of nucleons requires that the 
extraction of a gluon in the amplitude from a particular nucleon be matched by an identical gluon 
being extracted from the same nucleon in the complex conjugate. This along with the strong ordering 
of positions of the scatterings restricts the calculation to the set of nested scattering diagrams where 
the locations yi and are within the same nucleon. These diagrams are represented by the right panel 
of Fig. [5l Diagrams with more gluons per nucleon have smaller length enhancement factors, this is the 
reason that pomeron like exchanges which include at least four gluons per nucleon are suppressed in 
the large nucleus limit. Since the transverse size of a nucleon is much larger than the inverse transverse 
momentum carried by either exchanged gluon, one may integrate over the mean transverse location of 
the two insertions. This equates the transverse momentum that is exchanged between the jet and the 
nucleon, in the amplitude and the complex conjugate. For the z*^ scattering, this may be expressed as 



dy;dy';d%idY±{p\A'^{f±)A+{y'±)\p)e~''''^''^'''e'P^^^^^^ 



(27r)^5^(pi-p'^) / dyrdy'id'y^e 



-{p\A^{yj2)A^{-yj2)\p). (39) 



Given this pairing of interactions, the largest length enhancement that may be obtained is L" for n 
scatterings in the amplitude and complex conjugate. 

The other approximation is that the scatterings transfer small amounts of transverse momentum 
compared to the transverse momentum of the original quark, though both are at the same scale XQ. 
The possibility of scattering over n nucleons, greatly enhances the transverse momentum exchanged 
between the jet and the medium. Thus the transverse momenta from these scatterings cannot be 
ignored. However, since each is small compared to the transverse momentum of the jet, we may Taylor 
expand in them. In the case of spin independent cross sections, the leading term in the Taylor expansion 
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is the second derivative of the (5-function in terms of each of the small momenta k\ 



(40) 



The two factors of k]_ can be combined with the two gluon vector potentials in Eq. fl39l) to convert them 
into the correlator of field strengths {p\F~^^{y^)F^{y'j~)\p) . Higher transverse momentum derivatives 
will involve more factors of k± which may be converted through by-parts integration into derivatives of 
the field strength. Thus these are further suppressed without any extra length integration to enhance 
them. Thus we stop at just the double derivative in the Taylor expansion of the (5-function in terms of 
the transverse momentum. 

Incorporating these approximations, one obtains the differential hadronic tensor as, 



.(41) 



The latter quantity in square brackets is referred to as -D. This series may be resummed by noting that 
^ihi) ~ dW^'^ /cPlq. obeys the diffusion equation. 
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which has a normalized solution given as 
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(42) 
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(43) 



The transverse momentum squared gained by a parton as it traverses a length L is given as 
4:DL~ = qL~ . Thus the jet quenching transport coefficient q = AD. 



2.2.2 Multiple scattering of quark and radiated gluon 



With the derivation of the single gluon emission cross section and the multiple scattering cross 
section for a single quark without radiation, we can surmise the form of the multiple scattering single 
gluon emission cross section. Consider the general process depicted in Fig. El A quark is produced in 
the DIS on a large nucleus at location i/q. We consider the process where, in the amplitude, this quark 
will radiate a gluon without a re-scattering before the vertex at which the hard gluon is radiated. The 
produced quark and gluon will then scatter r and s times, respectively, on the incoming nucleons. The 
first location where the gluon scatters will be denoted as (c- The first location where the quark scatters 
will be denoted as yc There will be additional scatterings but the q = r + s scatterings are special as 
they represent the interference contributions. On the complex conjugate side, the produced quark will 
itself scatter r + s = q times on the same nucleons and then radiate the gluon. The location of the last 
such scattering which occurs on the parent quark in the complex conjugate and on either the produced 
quark or radiated gluon in the amplitude will be denoted as yE- Thus the radiated gluon is produced 
at two separate points: it is just after y^ in the amplitude and just after y'^ in the complex conjugate. 
After this cross scattering, the gluon will scatter m times and the quark n times in both the amplitude 
and the complex conjugate. The location of the first independent scattering on the gluon is at location 

and that on the quark at location yj- 
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Figure 6: DIS on a nucleus leading to the formation of a quark jet with a radiated gluon. 



The Feynman integral for this process may be expressed as (see Ref. [M] for a detailed derivation 
of this equation) , 
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(44) 



In the equation above q = r + s. While not explicitly mentioned, all scattering points on the same line 
(quark or gluon) are strongly ordered. The only unspecified orderings are scatterings on the gluon versus 
those on the quark after emission. This expression may be easily generalized from Eqs. (I26II25|38I) The 
(5-function 5{l_[_ + /gj^) is shifted to 5{lq_^J and is then broadened by multiple scattering as in Eq. (!38|) . 
There are now r + s + m + n two gluon matrix elements corresponding to the number of scatterings. 

While there is no simple means to derive the second line of Eq. (l44l) . it can be qualitatively understood 
as follows. Consider the case of gluon radiation without scattering as described in Sect. 12. 1[ The factor 
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of l/l^ in Eqs. fl25ll26p may be more naturally expressed as 

where is the transverse momentum of the gluon just after emission and e is the polarization vector 
of the produced gluon. This expression naturally arises in the case of one gluon bremsstrahlung from a 
hard quark by evaluating the square of the matrix element M. ~ J°''^A'^. This will be modified in the 
multiple scattering case. In A~ = gauge with the gluon moving in the (— )-direction, the dominant 
components of the projection vectors are the transverse components. In the case of no scattering the 
transverse momentum at emission is identical in the amplitude and complex conjugate. In the case 
of multiple scattering, the transverse momentum of the radiated gluon in the amplitude and complex 
conjugate is identical at the cut line but not at the point of emission. Thus, at the point of emission, 
the factor is modified as. 
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(46) 



Summing over the projections of the radiated gluons, yields the second line of Eq. (jSj). In the amplitude 
part (the right side of the expression), the gluon is radiated after the r + s scattering on the original 
quark. Of the transverse momenta gained by the parent quark in the r + s scattering, only a fraction 
y is transferred to the radiated gluon. This is the reason for the factor of y in the amplitude terms in 
Eq. (USD. 

The third, fourth and fifth line of Eq. (jUj) contain the products of the two-gluon matrix elements 
at the locations of the various scatterings. The first set on the third line contains the collection of 
correlated scatterings where the parent quark scatters in the amplitude at the locations i = 1 to g and 
the produced quark and gluon scatter off these same locations in the complex conjugate. The next set 
of scatterings on the fourth line at locations j = 1 to n represent the scatterings of the final produced 
quark in both amplitude and complex conjugate, i.e., there is no interference in these scatterings. The 
last set on the fifth line represent the independent scatterings on the final outgoing gluon line. 

The last two lines of Eq. (jUj) contain the phase factors that arise from the interference between 
the gluon emission from different locations in the amplitude and complex conjugate. The phase factors 
may be somewhat more straightforwardly motivated. There are in principle three types of vertices 
that will lead to the radiation of a hard gluon. These are illustrated in Fig. [71 The left-most diagram 
represents the case where the incoming quark is taken time-like off-shell, indicated by the filled dot. 
In the first term of the phase factors in Eq. (l44l) [sixth line of Eq. (144]) ] the quark in the complex 
conjugate is taken time-like off-shell by the scattering at y~^. This immediately introduces a phase 
factor of exp(— ip+XL?/^); the quark goes off-shell and immediately radiates a gluon. The 6'-function 
expresses the fact the the other propagators are all on shell and thus the location C,i is far ahead of ?/^. 

The second diagram in Fig. [7| expresses the possibility that the quark after the scattering remains on 
shell and radiates a space-like gluon which is then brought on-shell by the first independent scattering 
on the gluon. In the second term of the phase factor of Eq. (jH]), this occurs in the complex conjugate 
at location C,i. The exponential factor is exp(— ip+XL^/E). The negative sign is due to the fact that 
the propagator goes space-like. This in indicated by the hollow blob in the figure. Such a contribution 
occurs twice in the calculation: it occurs in the case where the initial quark propagator is taken on 
shell first or when the final quark propagator is taken on shell first. The fourth exponential factor is 
the instance of the latter pole in the complex conjugate. The third phase factor corresponds to the case 
where the on-shell quark radiates an on-shell gluon and goes off-shell. It is brought back on-shell by a 
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Figure 7: Diagrams which lead to LPM interference. 




Figure 8: Diagrams which are leading in the higher twist approach. 



hard scattering. In the case of the quark in the complex conjugate in Eq. (jM]), this occurs at location 
yj and has the corresponding exponential factor. The remaining phase factors for the gluon emission 
in the amplitude have similar explanations and thus similar phase factors. 

In the higher twist formalism, one takes the hard scattering hierarchy, in the sense that the hard 
scattering scale Q is assumed to be much larger than the medium scale ~ T or the scale of Debye 
screening. Thus, one tries to isolate terms which are minimally suppressed by the hard scale Q. This 
is obtained, as in the case of no emission, by Taylor expanding the small transverse momenta k\. The 
expansion of the 5-function can be shown to yield a similar diffusion equation as for the case of no 
emission. The emission yields a factor of the splitting function P{y) and phase factors which appear 
in the last line of Eq. ( 144|) . Of these, the phase factor which is dependent on the origin (e^^^^^^o ) 
represents the terms where the jet exits the hard interaction off-shell and immediately radiates a hard 
gluon. This is represented by the leftmost diagram in Fig. |8l Terms containing this factor in a product 
with the complex conjugates of similar diagrams as well as those where the hard parton is driven off 
shell later in its history are representative of vacuum like emission and vacuum medium interference. 
These are the leading contributions in the higher twist formalism. A subset of these are represented by 
the diagrams of Fig. |8l One also has to include contributions where the quark propagator after emission 
becomes space-like off-shell. 

Including these contributions along with the contributions where all three types of scattering may 
not be present and then Taylor expanding in we obtain the next-to- leading power suppressed 
contributions, i.e., those suppressed by one extra power of (or rather l\) as a sum of three terms 
each with a different power of the gluon momentum fraction y. The leading contribution to the radiated 
gluon spectrum in y arises from purely gluon scattering, and is given as 



We have dropped the subscript / from the position Cj in writing the above equation. The next con- 




(47) 
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tribution is proportional to y and includes the cross-scattering terms between the final radiated gluon 
and the quark before the split. This contribution to the radiated gluon spectrum is given as 



dydl\ 2-kI\ 



yP{y) dC'^[2-2cos{p+XLC)]- (4J 



The last contribution at this level of power suppression, suppressed by two powers of is the set of 
diagrams where the initial quark scatters in the amplitude and the complex conjugate. These diagrams 
were not included in the Eq. (l44l) or in the corresponding diagrams in Fig. [61 The contributions from 
these are 



dN^ a,C, 



fy'P{y) f dCCpq{C) [2 - 2 cos (p+x^C")] • (49) 
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Summing the above three equations yields the next-to-leading power correction due the single gluon 
emission cross section due to multiple scattering of the parent quark, radiated gluon as well as the 
produced quark. The argument for terminating at this order of power correction is the hard factorization 
limit of ^ qL~ , i.e., higher power corrections are further suppressed. 

2.2.3 Length dependence of energy loss and extensions to heavy-ion reactions 

In general, the application of Eqs. ( H7|48|49p to the computation of a medium modified fragmentation 
function tends to yield a complicated dependence on the length L~ of the medium. However, in 
earlier publications on the higher twist method, one presents a simple and approximate relation for the 
fractional energy lost in the small y limit to illustrate the length dependence of the energy loss. This 
has recently been criticized in Refs. [751 ES] • In the next few paragraphs we outline this estimate and 
review the critique against it. 

In the small y limit, one may focus solely on the first term of the differential gluon radiation spectrum, 
i.e., Eq. ( H71) . The fractional energy loss may then be estimated as. 



We say "estimated" as the integral over /_|_ is integrated from zero to infinity which is never carried out 
in realistic HT calculations. Note that both limits violate the basic assumptions of the higher twist 
method: Being an expansion in inverse powers of the hard scale l^, the radiated gluon's transverse 
momentum cannot be taken to zero without resumming all higher order contributions which have been 
dropped in arriving at Eq. f HTj) . Also the maximum radiated transverse momentum is limited by the 
virtuality Q of the jet, which in turn is limited by the energy q~ , i.e., l± < Q . Thus l± cannot be 
taken to infinity with q~ remaining finite. 

Ignoring these constraints, a simple minded estimate may be obtained by carrying out the l± inte- 
gration from to oo first and retaining only the leading term in y. We obtain, 

(,).r"*^#^r<-«r (51) 
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In a medium the range of y is limited by the fact that the momentum fraction xl = 2q-p+y{i-y) ^ 
i.e., a proton can at most impart its entire light-cone momentum to the jet in a single gluon exchange. 
Thus the y integration is not divergent. Assuming that q is independent of the location (~ we obtain 
the simple and well known relation that, 

{AE) oc qL~\ (52) 
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with the constant of proportionahty depending on the range of y integration. This, in turn, depends 
on the maximum of l\ = and on the energy q~ of the hard jet. Note that one cannot replace 
by oo in this part of the calculation. We should point out that this relation is quite approximate and 
in any realistic calculation the modification of the fragmentation function has a much more involved 
length dependence that what is suggested by the above equation. 

Hand waving estimates such as the above have come under some criticism recently for an entirely 
different reason. In obtaining the phase factors in Eq. (l44l) all factors of k± in the phases have been 
dropped. This is justified as kj_ dependent phase factors are generically of the form exp[i 2g-y(i-y) C~]; 
and differentiating these with respect to leads to multiplicative factors of the form C~/q~, which 
in the large q~ limit are small. This assumes that the coefficients of these terms are small. As the 
authors of Refs. [751 ES] have shown, this is no longer the case when one integrates over l±_ and allows 
l± to become large as in the estimate above. When l± is integrated up to oo the transverse momentum 
derivatives of the phase factors cancel those from the hard part. However in the limit that 2q'-'^{i-y) ~ ^ 
derivativatives of the phase factors are small and can be ignored. 

As pointed out previously, the HT approach is cast in the framework of DIS on a large nucleus. The 
medium in this case is confined. However the basic formalism may be straightforwardly extended to 
the case of a jet propagating through a deconfined medium. This extension is illustrated with the case 
of transverse broadening as in Eq. fHTj) . In this case the part of the hadronic tensor (or cross section) 
where the virtual photon strikes an incoming quark and sends it into the nucleus is entirely contained 
within Wq'^ . This part, which is specific to DIS, has been factorized from the remaining process which 
represents the multiple scattering of the hard quark in the medium. 

Extending this formalism to the case of jet propagation in a hot deconfined medium consists of 
replacing the hard cross section to produce a hard quark in DIS: 

EL^dao _ al^ L,^dWr _ aL Vl-^D 27rg^/.(^B )^^(V) ^^^^ 



d^L2d\jn^dy 2TTsQU%^dH^dy 2ttsQ^ d%^dH^dy 
with the cross section to produce a hard quark or gluon in a heavy- ion collision with a momentum pT, 
d'^a^ 1 /■ , ^ . I ^.daab^cxiQ,XaP,XbP,PT) 

- dXaGa{Xa,Q)Gb{Xb,Q) -p: . (54) 



dyd'^PT TT y dt 

What remains is the final state multiple scattering of the produced quark. This is assumed to be 
formally identical to the case of jet propagation in cold matter. The only change is that the expectation 
of the two gluon operator which is an input in the expression for the transport coefficient q Eq. ( HTll 
will now be evaluated in a hot deconfined medium, i.e., 



= %^ f 'JL±y±e-'[^'-'-'^) 5:(n|e-»F»-(,-.y,)F-(0,0)|n) (55) 

In the equation above, \n) represents a state in a thermal ensemble and "H is the hamiltonian operator 
and (3 is the inverse temperature. The exact value of q obtained depends on the methodology used in 
modeling the medium. For example, assuming that the medium can be represented as a weakly coupled 
plasma of quarks and gluons in the Hard Thermal Loop (HTL) approximation yields the expression in 
Eq. dZSI). 
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2.3 Other approaches and pictures of the medium 



2.3.1 AMY approach: Hard Thermal Loop field theory 

Unlike the HT approach, the formalism of jet modification based on the work of Arnold, Moore and 
Yaffe (AMY) starts with a rather precise definition of the medium. The medium is assumed to be 
composed of quark gluon quasi-particles with dispersion relations and interactions given by the HTL 
effective theory. Thus all quasi-particles in the medium have thermal masses ~ gT and their scattering 
is dominated by soft scattering. The hard jet is assumed to have a virtuality scale comparable to the 
Debye mass or thermal mass. Thus there is no jet-like hard scale in the problem. One constructs a 
consistent perturbation theory at the T — > oo limit and applies it to a realistic temperature. 

In spite of these differences, at the diagram by diagram level, the AMY approach evaluates diagrams 
which are topologically equivalent to those evaluated by the higher twist (HT) approach. To illustrate 
this we consider a generic diagram in the HT approach and convert it to a diagram similar to those 
evaluated in the AMY scheme. In the left top panel of Fig. [HI we consider the case that a hard quark 
produced in a hard scattering then radiates a hard gluon and both final partons scatter once off the 
medium. While in the case of HT, the dominant contribution contains the case where the initial quark 
could be quite off-shell and radiates immediately, this is not the case in AMY. The incoming parton is 
close to on-shell and will radiate a hard gluon only on being stimulated by scattering. As a result, the 
transverse momentum of the radiated gluon is always of the the order of gT. 




Figure 9: Starting with a higher twist diagram and obtaining the equivalent AMY diagram. 



As a result of the initial quark being on-shell, the momentum of the quark in the amplitude and the 
complex conjugate is identical, indicated by the double notches on these lines in the left top diagram of 
Fig. ini In this limit, the final state scattering part can be factorized from the initial structure function, 
as indicated by the right top diagram. The final state scattering part now contains an incoming quark 
which is on-shell and then is driven off shell by the scattering experienced by the remnant quark and 
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the radiated gluon. The complex conjugate also contains the same process. The blobs which connect 
to the scattering gluons can be replaced by the imaginary parts of the HTL self-energies. Given that 
the initial quark in the amplitude and the complex conjugate are identical, we can consider these to 
also arise from a cut diagram as shown in the bottom right panel. This final diagram corresponds to 
the type of self energy diagrams evaluated by the AMY scheme. 

While the diagrams considered by AMY and HT are topologically similar, the limit in which these 
are evaluated, as well as the pole structure, is quite different. As pointed out in the previous section, 
the leading contribution in the higher twist always contains, either in the amplitude or in the complex 
conjugate, a term which represents the quark being produced far off-shell in the initial hard scattering 
and immediately radiating a gluon. This is missing in the AMY approach. The dominant contributions 
in the AMY approach consist of interferences between different medium stimulated emissions along 
with the squares of the amplitudes themselves. As the single gluon emission cross section is evaluated 
in a different limit in the AMY approach, the formalism used to iterate this is also different and will be 
discussed in the next chapter. 

2.3.2 GLV approach: Opacity expansion 

This scheme named after Gyulassy, Levai and Vitev is very closely related to the HT approach. The 
primary difference is that the GLV scheme makes a very specific assumption regarding the structure of 
the medium. The medium is modeled as separated heavy static scattering centers with color screened 
Yukawa potentials. This model of the medium was originally introduced in the seminal work of Gyulassy 
and Wang [2U [251 126] and is often referred to as the Gyulassy- Wang model. The cross section for the 
interaction of the jet parton with one of these heavy medium partons is 

daei _ CrC2 

~ (27r)2(iV2 - 1) 

Where, Cr is the Casimir for the jet parton and C2 is for the medium parton. Alternatively, one may 
think of this as a specific model to evaluate the gluon field strength correlator of Eq. (ES])- The quantity 
/J, represents the Debye mass which screens the potential of the scattering centers and represents one of 
the parameters of the model. The other parameter is the mean free path of the jet A = (pcr)~^, where 
a can be obtained from integrating Eq. (!56|) and p, the density of scattering centers in the medium can 
be obtained from entropy considerations. 

Beyond this, the assumptions regarding the hierarchy of scales where the forward light-cone mo- 
mentum of the jet q~ ^ l± the transverse momentum of the radiated gluon [which in turn is much 
larger than the (— )-components] are the same as in HT. However, unlike the HT, there is no concept 
of scale evolution, thus the model of the medium as incoherent scattering centers is meant at the mo- 
mentum scale of the medium. The first term to be considered is the vacuum radiation term as in the 
HT scheme. This is followed by terms that involve the interference of medium induced radiation and 
vacuum radiation. As in the HT scheme, the leading contribution in the GLV approach also involves 
one scattering with the medium and is proportional to the number of scattering centers as seen by the 
jet. This is expressed as L/X, where L is the path length of the hard jet and A is the mean free path for 
a scattering. The former quantity depends on the medium in question while the latter is a parameter of 
the calculation. The ratio of L/X, called the opacity, is similar to the length enhancement factor in the 
HT approach and, combined with the cross section for a scattering, represents the expansion parameter 
for this calculation. 

Thus calculations which include only the vacuum radiation and one scattering, which are propor- 
tional to one power oi L/X and are called up to first order in opacity. The additive term in the cross 
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section to produce a single gluon at first order in opacity is 
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The major difference with the next to leading twist term in the HT expansion is the integral over the 
exchanged gluon transverse momentum k± and the explicit appearance of the the Yukawa cross section 
to scatter off a heavy scattering center. In the equation above, n[C,) = 6{L — Q/L is the longitudinal 
profile. 

Terms with two scatterings are proportional to (L/A)^ and so on. Note that leading order in opacity 
does not mean that the jet has scattered only once in the dense medium. Rather, it represents the 
case that at most one scattering was involved in a single gluon radiation. The basic methodology of 
the GLV approach is to compute the single gluon emission kernel, order by order in opacity. For dilute 
media or short path lengths one uses only the leading order in opacity, while denser or more extended 
media require higher orders. Higher orders in opacity have been explored analytically [IH], but it has, 
so far, not been possible to analytically resum the full series in opacity. 

2.3.3 ASW approach: Multiple-soft gluon exchange 

Similar to the GLV scheme, the ASW approach named after Armesto, Salgado and Wiedemann also 
assumes that the medium may be modeled as a collection of heavy static scattering centers with Debye 
screened potentials. In the calculation of the single gluon spectrum, the ASW scheme also invokes the 
small y limit where y is the longitudinal momentum fraction of the radiated gluon. As in the case of the 
HT scheme, it also finds that the leading contributions in this limit consist of diagrams where primarily 
the gluon scatters or diagrams with cross scatterings where the quark in the amplitude scatters and 
interferes with gluon scattering in the complex conjugate. 

In the ASW approach, the multiple scattering analysis is carried out in impact parameter space. The 
square of the amplitude where the gluon scatters m times in both amplitude and complex conjugate, is 
at a transverse position in the amplitude and at a transverse position y± in the complex conjugate 
is given as 



where, Q and C,f indicate the longitudinal positions over which the scatterings occur. The quantity n{Q 
is the opacity at the location C, and a is the dipole cross section of two gluons separated by a transverse 
distance u± — y±. The gluon propagation from Q to in the amplitude represents one part of the 
dipole, where as the propagation of the gluon in the complex conjugate is the other part of the dipole. 
The case of multiple scattering is obtained by summing over m to obtain the obvious exponential with 
the argument in parenthesis in Eq. 0581] . 

In the region before (j there is cross scattering with scattering on the incoming quark in the amplitude 
matched with scattering on the gluon in the complex conjugate and vice versa. This case is more 
complicated; starting at the location (j, the quark propagates in the amplitude and a quark gluon 
pair in complex conjugate. The result is given by a path integral where the dipole moves in the two 
dimensional space of relative transverse location. The diagrams involved in this part are near identical 
to those of the AMY scheme. However, in the AMY case, the entire calculation is carried out in 
momentum result all positions are integrated from — oo to oo. As the calculation here is 

carried out in position space, the numerator factor of lj_ ■ in Eq. (l46l) is replaced by the appropriate 
derivatives over transverse position. The final expression for the cross section to radiate a gluon with 





30 



transverse momentum l± and energy u and momentum fraction y (or alternatively the distribution of 
radiated gluons in transverse momentum and energy), is given as 
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dxdt^i 
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For a short range potential one may take the saddle point approximation around the — limit by 
replacing, 



(60) 



This approximation allows the above path integral to be carried out. While introduced in this way, 
the one transport coefficient q has the standard meaning as in any other formalism: the transverse 
momentum gained by a hard parton (in this gluon) per unit length. 

Unlike the HT or the GLV approach, there is no explicit introduction of the vacuum contribution 
or vacuum medium interference. The fact that a hard jet in DIS or a heavy-ion collision is not created 
at yl^ = — oo is incorporated by introducing a finite starting point to the longitudinal integrals. In 
this case starting yf integral at Co" usually set to be the origin. There is no systematic insistence that 
the off-shellness from the vacuum contribution be much larger than that from the medium scatterings. 
Given all these approximations, the single gluon cross section from multiple scattering is found to be 

dependent on only one parameter: Wc = f d((q{() = -qL"^, as (in the UcL — > oo limit). 
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2.4 Parametrizing the effect of the medium: dN/dy, T and ag 

All the different schemes above depend on one or two parameters. In the case of the HT, the unknown 
quantity is the two gluon field strength correlation in the medium in question: 

I dy~d^y^e-'^'-^'--'-{F^^{y-,y^)F;m. (62) 

Expressing F^^^ = F'^^^Vj^ where Vy is the light-like velocity of the hard parton we note that the above 
quantity is simply the Lorentz force correlator in the system. Since this quantifies the transverse 
force experienced by a jet as it passes through a medium, it may be straightforwardly related to the 
transport coefficient q as in Eq. (!55|) . Since k_[_ ^ ^qcd, the integrals over transverse location y±_ are 
rather limited. On the other hand = k\/2q~ may well be non-perturbative, so that the y~ integral 
is limited by non-perturbative effects, such as the the confinement length in cold nuclear matter and 
the screening length in hot deconfined matter. 

To obtain some form for the two-gluon field strength correlator in Eq. fl62|) . one requires a microscopic 
model of the medium. The situation in the ASW approach is similar. Although a specific medium model 
of static scattering centers is used in the derivation, the final result depends only on the transport 
coefficient g, suggesting that the energy loss within the multiple soft scattering approximation (i.e. 
when the saddle point approximation of Eq. f l^U]) is valid) does not depend on the specific details of the 
medium, but measures directly the mean transverse momentum exchange per unit path length. 

In the AMY scheme, which is cast completely within the HTL formalism of finite temperature field 
theory, the scattering is not factorized from the gluon emission process. Since HTL perturbation theory 
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is cast in the limit that there exist separate scales T, gT, the two obvious parameters of the theory are 
T and g (or rather a^). Note that this formalism may only be applied to hot deconfined matter. The 
coupling constant is the same for the jet and for the medium. It is allowed to vary as a fit parameter 
and set by fitting one data point. 

Since the AMY formalism makes a precise assumption regarding the microscopic structure of the 
produced matter, it provides a well-defined way to calculate the the transport coefficient q, given a T 
and an Og- Following Ref. [77] we note that q may be defined as, 

where, dT/(Pk± is the differential rate for elastic collisions of a hard parton with quasi-particles from 
the medium. While the expression for the differential rate has slightly different forms depending on 
whether the exchanged gluon is soft or hard, we will approximate the rate with the expression for soft 
momentum transfer: 

jfT Cn g'Tml 

where mz) is the Debye mass, m|, = ^^"g"^^ iyNc + Substituting the expression for the differential 

rate in Eq. ( l63l) . one may easily calculate q given the range of integration of the exchanged momentum. 

One notes immediately that the expression for q is bounded from below and thus the lower limit of 
integration may be taken to zero. This is often done in many different approaches. The upper limit is 
bounded by kinematics of the process, i.e, the outgoing two partons may not carry more energy than 
the incoming partons. This constrains the upper limit to q^^^^ ~ \/ AlET where E is the energy of the 
hard parton. The energy loss calculation in the AMY is based on rate equations similar to Eq. f lM|) and 
thus uses the same transport properties, but without explicitly calculating the transport coefficient q. 

In the case of the GLV approach, one needs to set two parameters: the Debye screening mass and 
the mean free path of the radiated gluon A. The mean free path may be obtained using A = l/(pcr), 
where p is the density of scattering centers and a is the integrated cross section for a gluon to scatter 
off these. Because of the dependence on p, the total number of gluons at mid-rapidity dN/dy has often 
been quoted as a measure of the medium density in GLV calculations [75]. However, it should be noted 
that the Debye screening mass and thus the temperature also enters in the GLV formalism via the 
scattering potential Eq. fl56|l and the cross section a. 

In conclusion, this means that in all four formalisms, a microscopic model of the medium is needed 
to set the input parameters for the energy loss calculation, be it the two gluon field strength correlator, 
the transport coefficient g, or m^) and A. The AMY approach makes a specific choice for a HTL plasma. 
Given such a model, the local medium parameters can be related to single quantity in the medium, 
such as the temperature T in HTL. It is, however, worth to keep in mind that the radiative energy 
loss in all models is driven by transverse momentum exchanges, so that the energy loss is governed by 
the transport coefficient g or a closely related quantity. In addition to the microscopic model of the 
medium, a realistic calculation of energy loss in heavy collisions also requires a macroscopic model of 
the medium, specifying the space-time dependence of the local properties of the plasma. This aspect 
will be further discussed in Chapter [31 
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3 Incorporating Multiple radiations 



In the preceding chapter we described the calculation of medium induced single gluon radiation from a 
hard parton traversing a dense medium. Imagine a gluon radiated with a transverse momentum and 
a forward momentum l~, its formation time is given as, r = l~ For realistic values of = 3 GeV^ 
and /~ = 6 GeV say from a jet with forward momentum q~ ~ 20 GeV, we obtain a formation time of 
r ~ 0.4 fm. Thus in almost any medium with a length L > 1/2 fm there will be more than one such 
radiation. In this chapter we describe the different ways in which the single gluon emission is iterated 
in different jet modification schemes. 

The incorporation of multiple emissions completes the theoretical calculation of radiative energy 
loss for light partons. The final form of the result depends on the observable in question. For the com- 
putation of the single hadron inclusive cross section one computes the medium modified fragmentation 
function D^{z), where z is the momentum fraction of the hadron h with respect to the momentum 
of the originating hard parton. We should point out that, unlike the case in vacuum where the hard 
parton produced in a hard collision continues to monotonously lose virtuality till the onset of the 
non-perturbative process of hadronization, the process of multiple radiation in a medium may be far 
more complicated, including a series of both virtuality increasing collisions and virtuality decreasing 
emissions. To date, the space-time structure of this process has not been completely elucidated. The 
different approaches of the various schemes represent approximations to the actual mechanism. 



3.1 Higher twist: In medium DGLAP 

A hard parton from a large collision is considerably off-shell when it is produced. This parton 
then proceeds to lose this virtuality by a series of gluon emissions ordered in virtuality or transverse 
momentum. The radiated gluons are themselves virtual and continue to radiate lower virtuality gluons. 
The progress of this shower may be described using pQCD as long as the virtuality n is large compared to 
Aqcd- Beyond this, one describes the process with a phenomenologically fitted fragmentation function 
D^lzjfi"^) which contains all processes up to the scale /i^ which produce at least one hadron with 
momentum fraction z. The effect of multiple gluon emissions from the scale up to some predefined 
virtuality or scale is incorporated by evolving the fragmentation function using the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation from /i^ up to Q^. 

The DGLAP evolution equation may be understood as the resummation of an arbitrary number of 
ordered emissions which connect a hard parton with a virtuality with a parton with lower virtuality 
/i^. For the non-singlet fragmentation function D^^s = ~ D^i which is always vanishing when the 
fragmenting parton is a gluon, this process may be straightforwardly expressed as. 



Diz,Q') = D{z,f^') + ^ dli^:i^^D[-,f,'] (65) 



We have ignored the ordering in the momentum fraction y and have only written the first three terms in 
a series where each term contains a product of integrals over transverse momenta where the transverse 
momenta at each emission is limited by the preceding radiation. The factor P{y) represents the q qg 
splitting function, i.e., the probability for a quark to radiate a gluon and retain a fraction y of its 
momentum. This series can be resummed by noting that differentiating with respect to log(Q^) yields 
an identical series as above for the shifted fragmentation function D{z/y), convoluted with the splitting 
function P{y): 

dD{z,Q^) 



dlogQ^ 



(66) 
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For an arbitrary fragmentation function, the evolution equation for the quark-to-hadron fragmentation 
function is coupled with the evolution equation for the gluon-to-hadron fragmentation function. 

The jet modification formalism which is closest to this line of argument is the HT formalism. It 
assumes that while scattering processes may raise the intermediate virtuality in a given amplitude prior 
to or just after a gluon emission, there is an overall drop in the average virtuality between the incoming 
parton and the outgoing pair of radiated gluon and remnant parton in the single gluon emission cross 
section. Note that multiple scattering tends to progressively raise the virtuality of the propagating 
partons. Thus the above assumption is valid in the regime of densities and incoming virtualities where 
the formation time of the gluon is so short that the drop in virtuality from the emission dominates over 
the rise due to multiple scattering over the same time, i.e., 

gr = ^</i^/i>V^. (67) 

We also note that gluons with a formation time r longer than the length of the medium L are not 
influenced by energy loss, which sets a minimum virtuality: t < L ^ l\ > 21' /L. Thus for the 
formalism to remain consistent till the exit of the jet from the medium we require the minimal condition, 

VW =^^AE = -qL\ (68) 

The last equality is obtained by relating the light-cone momentum of the radiated gluon with the 
light-cone loss and thus the energy loss. 

In a calculation one simply replaces the vacuum splitting in Eq. ( 166|) with the in-medium splitting 
function derived in the section above. It is assumed that the hard parton exits the medium with a 
virtuality (related to the transverse momentum of its radiations) //^ ~ 2Ey/ L ^ EjL^ where E is the 
energy of the initial hard parton, L is the length of the medium and y is a representative momentum 
fraction carried by the radiated gluon; this is approximated as y ~ 1/2 in the second part of the equality. 
As pointed out earlier, a dynamical calculation of the loss of virtuality with emission as a function of 
the distance travelled by a hard parton is still lacking. 

Given the measured fragmentation function at the scale /i^, this is evolved up to the hard scale of 
the process which represents the largest possible initial virtuality. Results of such an evolution in a hot 
QGP-brick will be presented at the end of this chapter. The result of the evolution is not very sensitive 
to the choice of the upper limit of the evolution; for instance, using (5^/4 will yield very similar results 
as long as ^ n^. 

3.2 ASW and GLV: Poisson convolution 

The ASW and the GLV schemes share a common microscopic picture of the medium: that of heavy 
static scattering centers. The also use an almost identical methodology for iterating the single gluon 
emission. Unlike the HT which computes the change in the distribution of hadrons (or the AMY scheme 
which computes the change in the distribution of partons), the central quantity that is calculated in 
these formalisms is the probability distribution P{AE) for a hard parton to lose an energy AE, via an 
arbitrary number of gluon emissions. This probability is then used to shift the fragmentation momentum 
fraction z and define a medium modified fragmentation function. 

Contrary to the HT approach, in these formalisms, the virtuality evolution of the parton is not taken 
into account. A hard jet propagating through the medium will progressively gain virtuality by collision 
and lose it by emission. Each of these emission processes makes a minimal shift in the energy of the jet 
and to leading order one may ignore this shift, this is often referred to as the eikonal limit. Thus all 
these emissions may be considered as independent of each other. The fate of these emitted gluons is 
also not considered, the focus is solely on the emitting hard parton. 
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In either formalism, one calculates the differential distributions of radiated gluons, e.g., in Eq. (159|) . 
To obtain the differential energy distribution of radiated gluons, these expressions have to be inte- 
grated over the transverse momentum l±. The formalism does not have a natural cut-off for trans- 
verse momentum integration, but instead uses kinematic considerations to limit these integrations, e.g., 
^Max <]2iin{|/, 1 — y}E as used by the GLV scheme or in the small y limit l^""^ < yE as used by 
the ASW scheme. The different choices of kinematic bounds, along with slightly different definitions 
of y (light-cone momentum fraction in GLV and energy fraction in ASW) lead to somewhat different 
results [79] . 

Having limited the transverse momentum integrations, the differential energy or momentum fraction 
distribution dl /duj,dl/dx is well defined. The integral over this quantity yields the mean number of 
gluons radiated from the jet: 

.r f , dl f dl 

N„ = / dx— = / du—. 69 
J dx J du 

Given that each of these emissions are independent, the number of radiated gluons n is assumed to 
follow a Poisson distribution, giving the following probability distribution for energy loss: 
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Using the expression for dl /du the above probability distribution may be numerically calculated. Having 
lost energy A£', the jet now fragments a hadrons with momentum ph off a hard parton with energy 
E — AE. The shifted or medium modified fragmentation function is defined as. 



D{z,Q') = dAEP{AE)^—^^j^. (71) 

Under the assumption that the virtuality gain from the medium is comparable to that brought in by 
the hard parton, the scale used in the final fragmentation function is the same as the hard scale of the 
process. 



3.3 AMY: rate equations. 

In the AMY formalism, the hard jet is treated similar to a hard parton in the medium. The virtuality 
is assumed to be the same as that for any hard thermal parton in the medium ~ gT. In this sense it 
is similar to the assumption made in the ASW and GLV calculations. In this formalism, the radiation 
of a gluon by a quark or a gluon is an order as suppressed process and thus it occurs rarely enough 
that one may consider sequential emissions to be independent of each other. The rates for a quark to 
decay into a quark and gluon, a gluon to decay into a quark anti-quark pair etc., are now used to set 
up a Fokker-Planck equation which describes the change in the distribution of hard partons with time. 
There are two sets of equations which describe the change of the distribution of the sum of quarks with 
anti-quarks and gluons. For example the q + q distribution may be described as 

^ ^ ^ - . .... . (r.) 

In the equation above, the first term on the r.h.s. describes the process where a hard quark (antiquark) 
with a momentum p+l decays in to a quark (anti-quark) with momentum p and a gluon with momentum 
/. Hence, its contribution to the Fokker-Planck equation is proportional to the distribution oi q + q 
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at a momentum p + I. The second term on the r.h.s. represents the decay of the quark (antiquark) 
with momentum p into lower momentum quarks and gluons and thus represents a depletion of the 
distribution at that momentum. The last term proportional to the population of hard gluons with 
momentum p + / represents the process where the gluon decays into a quark anti-quark pair where one 
has a momentum p and the other has a momentum /. The factor of 2 accounts for the case where the 
momenta of the two fermions are reversed. 

The integrals over / run from — oo to oo. When the momentum of the outgoing gluon or quark (/ 
or p — Z) becomes soft, of the order of T, they will encounter Bose enhancement or Pauli blocking from 
the medium. These factors of distribution are already included in the expression for the rate dT/dldt. 
Negative values of I indicate absorption of thermal gluons from the medium, the distributions for these 
arise simply from the sign change in the functions which lead to Pauli blocking or Bose enhancement. 
The AMY formalism is different from the other three approaches as it naturally incorporates feedback 
from the medium which is missing from all other approaches. A similar, but slightly more complicated 
equation may be written down for the gluon distribution [SU]. These two equations are solved in tandem. 

In this set-up one represents an incoming quark jet by setting the distribution at t = to be a 
(5-function, Pq^q{p,t = 0) = 6{p — Ejet)- The gluon distribution at this time is set to zero. These 
distributions are then evolved up to the time of exit by means of the coupled Fokker-Planck equations. 
In spite of the completeness of this approach, it does suffer from the rates being computed at a fixed 
temperature. Which implies that the medium remains static during the formation time of the radiated 
gluons. As a result, the application of this formalism to short or rapidly varying media is somewhat 
tenuous. 

The final medium modified fragmentation function is obtained by convoluting the final distribution of 
hard quarks, anti-quarks and gluons with the respective shifted fragmentation functions, using Eq. fl7T|) 
with the probability of energy loss replaced by the distribution of hard partons. In most applications 
so far the scale at which the fragmentation functions are evaluated is chosen to be the hard scale of the 
process. 

3.4 Comparison of the different energy loss formalisms 

In the previous sections, we have outlined the calculation of the medium modified fragmentation func- 
tions within the four different schemes that are currently in use. We now attempt a standard comparison 
between these different formalisms in an identical medium. In a recent effort by the TECHQM collab- 
oration [81] this problem has been formulated as a computation of the modification of a single quark 
jet propagating through a homogeneous static medium of fixed length held at a constant temperature. 
This medium is often referred to as the "QGP brick", as temperatures are usually chosen to be high 
enough that deconfinement has set in. 

In order to carry out computations in all four formalisms, their medium parameters have be defined 
in this standardized medium. We present our own particular prescription of a standardized medium 
in this review. One starts out by setting the temperature T and a value of the in-medium coupling 
as = 0.33. Given these two parameters, calculations of jet modification may be carried out in the AMY 
formalism. Given these two parameters, q can be computed using Eq. fl63|) . Carrying out calculations 
for a quark-gluon plasma with 3 flavors of quarks {Nf = 3) using the HTL approximation for small k±, 
Eq. (IMl) . we obtain. 



The value of q obtained by the prescription above for the AMY scheme may now be directly sub- 
stituted in the HT and the ASW approach. This will in some sense place these three calculations on 




(73) 



with qf^^ = Vet. 
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the same footing, by using the q for a weakly coupled plasma of quark gluon quasiparticles in all three 
formalisms. 

The GLV formalism requires the specification of the Debye mass m£, and the mean free path A. The 
Debye mass is taken directly from HTL field theory m|, = i^Nc + ^ j . The mean free path A is 

calculated from the scattering rate 

A - J ^ ^''^ 

using the same HTL approximation for small k±, Eq. f l6^ . Unlike the calculation for q, the equation 
above has an infra-red divergence. This is controlled by restricting the lower limit to m/j, with the 
result: 

1 r"^^'' dr^ 

with the further approximation g^^^ » nir,. 

Figure fTOl shows a comparison of quark fragmentation function ratios using four different formalisms 
to calculate energy loss (HT, AMY, WHDG radiative and ASW) for a uniform medium with 2 different 
path lengths L and temperatures T in the regime relevant to RHIC. The Higher Twist (HT) uses q and 
L as input to calculate a fragmentation function in the medium. The Figure shows the ratio between 
the medium-modified fragmentation function and the fragmentation function in the vacuum. 

In the AMY formalism, the Fokker-Planck equations are solved to calculate the distribution of 
outgoing quarks and gluons starting from a mono-energetic distribution of quarks, the incoming quark 
with = 20 GeV. The outgoing quarks are then convoluted with the KKP fragmentation function to 
give a medium-modified fragmentation function. Gluon fragments were not included in this calculation 
for the purpose of comparison to the opacity expansion and multiple soft scattering approximations. 

For the multiple soft scattering approach, ASW, the publicly available quenching weights code was 
used to calculate the energy loss probability distribution P(A£') given the transport coefficient q and 
medium length L. The medium modified fragmentation function is calculated by convoluting the energy 
loss probability distribution with KKP fragmentation functions using Eq. (ITTjl . Note that only quark 
fragments are taken into account, fragmentation of the radiated gluons is not calculated. 

For the opacity expansion, we used the DGLV expressions from Appendix B of Ref. [78], with an 
in- medium gluon mass of nig = mo/\/^ and a quark mass niq = m£)/2. The Debye screening mass 
and 1/A are defined above. Multi-gluon emission is treated using a Poisson convolution for multi-gluon 
emission. The medium modified fragmentation function was calculated using Eq. fl7T]) . 

Comparing the fragmentation functions in Fig. [TOl we see a characteristic rise in the ratio at low 
z for the HT calculation that can be identified with gluon fragments (including secondary medium- 
induced radiations), which are not included in the other calculations. Because of the steeply falling 
parton spectrum, the large z behaviour is most important for high-p^ hadron production. For short 
path length L = 2 fm, HT, GLV and ASW give very similar results at high z > 0.5, for both T = 250 
MeV and T = 350 MeV. The AMY calculation shows a much larger suppression. This may be due to 
the absence of vacuum-medium interference in AMY. 

At larger L = 5 fm (bottom panels of Fig. [TOl) . larger differences between the various formalism 
are seen. The magnitude of the suppression in the GLV opacity expansion and the ASW multiple-soft 
gluon emission formalisms are similar, but the z dependence is different between the two. The HT 
result approaches the AMY result for L = 5 fm and T = 350 MeV. It should be noted, however, that 
in this case, the suppression becomes very large, more than a factor 10 at 2; > 0.5, thus violating the 
AE « E limit in which all the formalisms have been derived. 
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Figure 10: Comparison of quark fragmentation function ratios using four different formalisms for a 
uniform medium with L = 2 fm (upper panels) and L = 5 fm (lower panels). For both upper and lower 
panels the left plot is at T = 250 MeV and the right plot is at T = 350 MeV. For details, see text. 
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4 Heavy flavors and elastic energy loss 



In the preceding chapters we exphcitly discussed the case of hght quark and gluon energy loss. It is 
now almost established that for light partons, transverse scattering induced radiative energy loss, plays 
a dominant role (see Section l6.4.5p . There is, however, an entirely orthogonal mechanism: elastic energy 
loss. While the nomenclature is somewhat confusing, this refers to a process where the jet transferred 
energy and momentum to a constituent of the medium and the interaction with the medium did not 
stimulate the jet to radiate a hard gluon, i.e., it was not inelastic for the jet. On the medium side, the 
effect may be both elastic or inelastic. We should point out that existing calculations of elastic energy 
loss [H21 ESI [Ell [39] have focused specifically on 2-to-2 scattering, where the medium parton also stays 
more or less on shell. 

For heavy quarks, radiative energy loss is much reduced, due to the dead cone effect, which suppresses 
coUinear radiation. As a result, elastic energy loss may be more important for heavy quarks than for 
light quarks. In the remaining part of this chapter we will review the dead cone effect and how it 
reduces the energy lost by a hard heavy quark. Following this we will consider the inclusion of elastic 
energy loss within the jet modification formalism. 



4.1 Dead cone effect in radiative energy loss 



Imagine a hard heavy quark produced in a DIS event about to radiate a hard gluon. Physically this 
will occur in a process where a hard heavy quark in the nucleon (produced in large fluctuation) is 
struck by the virtual photon. An alternative to this process is when the incoming electron encounters 
a charge current interaction via a W~ exchange which converts a u quark in the nucleon to a 6 quark. 
In the interest of simplicity and given the factorization of final state effects from the hard cross section 
and the initial PDF, we ignore the specifics of the charge current interaction and only focus on the final 
state radiation of a hard gluon from the produced heavy quark. 

The simplest means to see the effect of the mass of the heavy quark is to introduce a mass M in the 
fermion propagators in Eq. fl2^ . This converts the equation to 
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Carrying out the integrations over z and z' yields the usual constraints, Ig = po + q — I and po = p'q- 
These are used to perform the p^ and Iq integrals. The difference with the case of the light quarks are 
essentially the factors of M that appear in the denominators, in the Dirac matrix structure and in the 
on-shell constraint 6{lg — M^). The argument of the 5-function may be simplified as (using Ig = p^ + q — l 
and / is on-shell), 

2p+q- + 2q+q- + 2p+r + 2g+/~ + 2g~/+ - = 

2P+q~ (^Xo(l -y)- XBil - y) - ^ - M^^ = 2P+g-(l - y) [xq -Xb-Xl- Xm] ■ (76) 
Where, xl = l\/2P^q~y{l — y) has the usual definition, while, xm = M^/2g^P+(l — y). 
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To simplify the trace, we take the small y limit i.e., the soft gluon limit, which reduces this expression 
to a form similar to that of Eq. f[T^ . In this limit, the only projection of the gluon propagator that has 
to be retained is G^^ = l\/{q^_)- As a result, one only needs to retain the terms proportional to 7" in 
the gluon emission vertices. This simplifies the trace over Dirac matrices: 



Tr 



(77) 



With these simplifications, the operator O^^ [defined for the case of light quarks in Eq. (123|) ] which is 
directly proportional to the cross section may be expressed as. 



000 ^ ^ dyt-—ll—-. (7J 



-.00 _ ( ^sCf f ,2 

27r J y{ll + ym^f' 

Defining the angle of radiation as 6' = l±/uj~ where u~ = q^y, and the mass dependent angle M/u~ as 
9o we obtain the gluon radiation distribution to be 



1 + 7^ > (79) 
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where dNg^gg/dO"^ is the angular distribution from a light quark which has a logarithmic divergence as 
6* —7- 0. The angular distribution of gluon radiation from the heavy quark is cut off at small angles. This 
is called the dead cone effect [SB] • The origin of this effect is essentially the shielding of the coUinear 
singularity by the large mass of the quark. Since gluon radiation from quarks occurs dominantly at 
smaller values of l±, the suppression of radiation in this region by the mass of the quark, leads to a 
smaller amount of radiative energy loss from a heavy quark compared with a light quark. 

In the case of gluon radiation from multiple scattering there will be multiple instances of the heavy 
quark mass as there will be multiple quark propagators. So far, the calculation of medium induced gluon 
radiation from a heavy quark has only been performed up to one scattering per radiation in the HT 
formalism [86l [87] , up to 1*^* order in opacity in a variant of the GLV approach (called the DGLV) [50] 
and in the multiple soft scattering approximation of the ASW approach ^5\. As in the case of light 
flavors we will follow the derivation in the higher twist approach. The expression for the radiated gluon 
distribution from a massive quark is given, following Ref. [86] as a shift in the formation time r/ of the 
radiated gluon and an overall multiplicative factor / as, 

'7"f = -p^ ^TTTT; and / = J. (80) 

^ ll + y^M^ [11 + y^M^f ^ ^ 

In both the GLV and the HT calculations, the above factors lead to a considerable reduction in the 
amount of energy loss by radiation of gluons. This has led to the realization that other sources of energy 
loss, such as elastic energy loss, which is a smaller effect for light partons, may be more important in 
the case of heavy quarks. 

While the radiative energy loss calculations in the GLV, HT and ASW schemes have all been modified 
to incorporate the effect of the heavy quark's mass, it was found that in all three cases, purely radiative 
processes cannot describe the observed suppression in the yield of non-photonic electrons from D and 
B meson decay (see also Section r6.4.7p . Attempts to incorporate elastic energy loss within the energy 
loss formalism and describe the heavy flavor suppression have been carried out by the DGLV and the 
HT schemes. In the last section of this chapter we will describe the inclusion of elastic energy loss in 
jet modification calculations. 
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4.2 Longitudinal drag and diffusion: other transport coefficients 



In the derivation of the multiple scattering of a hard parton and the ensuing gluon radiation in 
Chapter 2, we specialized to the case where the negative light-cone momentum of the exchanged gluon 
was assumed to be vanishingly small, i.e., k~ — )■ 0. All factors of the exchanged negative light-cone 
momenta k~ were ignored from the expressions and this led to the restriction of the entire process to the 
y'^ = plane. In reality, depending on the structure of the dense medium with which the jet interacts, 
one can have a non-vanishing amount of k~ exchanged. This is true not only for a massive quark but 
also for light partons. In this section we derive the leading correction to the differential cross section for 
a hard light quark traveling through a large nucleus encountering multiple scattering which may change 
both its transverse as well as its longitudinal momentum. These may then be easily extended to the 
case of a heavy quark. The final expression will be written in a factorized form where the hard dynamics 
of the quark propagation will be separated from the in-medium expectation of gluonic operators which 
will codify the respective in-medium transport coefficients. Once defined in such a factorized form, the 
results may be immediately extended to the case of jets in a heavy ion collision by replacing the hard 
production cross section in DIS with that in A-A and replacing the transport coefficients with those 
calculated in a hot deconfined environment. In real comparisons to date, these are often fitted to light 
hadron production data points. 

The reader may wonder about the focus on light cone loss as opposed to actual energy loss. This 
emanated from our insistence on defining as boost invariant a quantity as possible. For a large magnitude 
of energy loss in an extended time At, the measure AE/At changes with boost, but the light cone ratio 
is boost invariant: 

Aq~ _AE- Ap, 7(1 - (3)iAE - Ap,) 



AL- At-Az -f{l - P){At - Az) ' ' ' 

Starting again from the case of DIS on a large nucleus, we now define the triple differential distri- 
bution for the hadronic tensor for the case of no rescattering and no final state radiation where a quark 
with xb momentum fraction is struct by a hard virtual photon and sent back through the nucleus as 

= WrS%m-, - q-). (82) 

where W^^ = —g^'^'^.n^^Q'^^f^^XB)- The two (5-functions indicate that the three momentum compo- 
nents of the outgoing quark have a narrow distribution around /~ ^ q~, the large light cone momentum 
of the incoming photon, and lqj_ ~ 0. Being close to onshell, ~ l\/2l~ 

Considering the case of one scattering, the modified differential distribution may be expressed as 

dl~dn^ , _ , + ,2 dk~d^ki_{2Tif5{l~ -q- - k-)6\t± - ki_) TrM 



W^^= Wrg'J ^^^dY~dy-dyU'y^ 



Tr 



X 



7V <! (r+ r ) + 7T7^^^- f^± \ 1^ 



(27r)3 2{q- + fc") A^^ 

k"^ 

exp -i-^{y~) + iy± ■ k± - iy-^k' 
L 2g- 



X (P|A^(F-+y)A^(F-)|P). (83) 

One immediately notices the shifts in the three (5-functions which describe the final outgoing quark 
distribution. These are, as expected, proportional to the expectation of the two gluon operator in 
the nucleus. The Dirac structure is obvious from the fact that the quark immediately after the hard 
scattering on the virtual photon is almost on shell in both the amplitude and the complex conjugate 
and has a magnitude only in the negative light cone direction / = (0,g",0,0). The cut line emanates 
after the scattering and has a momentum [k~^,q~ + k~,k±\ with k~^ = k']_/2{q~ + k~). We make the 
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same set of approximations on the gauge fields as for the case of transverse broadening, i.e., A'^ ^ A± 
(recall that we are calculating in A~ = gauge) as in Section 12.21 

Given the scaling of the A""^ field, we obtain that the dominant contribution in the Dirac trace as 
given by 

Tr [7+7"7+(g~ - A;")7"] • (84) 

In the case of transverse broadening, one Taylor expands in k±; to obtain the distribution in light 
cone momentum one Taylor expands in k~ . In the interest of simplicity, we will ignore the factors of 
transverse momentum and simply focus on the light cone components. Taylor expanding the hard part, 
we obtain that the leading terms emanate from the expansion of the ^-function, 

S{l--q--k-) = S{1- - q-) - ^^^^^ ~ '^"^ (fc-) + \-^Al- - + • • • • (85) 

In both the current case of longitudinal drag and diffusion, as well as in the case of transverse broadening 
in Eq. ( I40l) . the 5-function is meant in the limit of a narrow Gaussian distribution and should not be 
thought of as a singular distribution. The small shifts k~ are infinitesimal compared to the jet energy 
q~ . In the limit of multiple such interactions we, once again, appeal to the central limit theorem and 
compute only the mean and the variance of the distribution. 

Unlike the case of transverse broadening there is no cylindrical symmetry to insist that the term 
proportional to k~ is vanishing. One thus obtains two sets of terms involving derivatives of the 5- 
function, one which involves single derivatives, which lead to a drag effect on the light-cone momentum 
distribution and the other which depend on double derivatives which lead to a diffusion in the light 
cone momentum. The coefficients of these terms will yield the elastic light-cone loss per unit light cone 
length e and the diffusion in light-cone loss 62. 

The coefficient of the the single derivative of the (5-function may be simplified as, 

C, = j rfy-d\^^^-^e-4^"-^'="^"+^'=-^-^(P| [^^-A^^{y)A'^^{Q)] \P). (86) 

In the equation above \P) represents the nuclear state. We have dropped the Y~ in the arguments of 
the gluon field based on translation invariance in a large nucleus. Color confinement insists that the 
quark struck by the hard virtual photon emanate from the same nucleon in both the amplitude and 
complex conjugate. Similarly the gluon off which the outgoing quark scatters must also be restricted to 
the same nucleon in both the amplitude and complex conjugate. In the model of a nucleus made up of 
almost free nucleons (in the high energy limit), the correlation between these two nucleons is p{y~)/2p'^ 
where p is the nucleon density evaluated along the path of the propagating out-going quark and p'^ is 
the mean light cone momentum of a nucleon in the nucleus = P~^/A. 

Incorporating the above simplifications one obtains the light-cone loss per unit light cone path 
travelled by a struck quark without radiation in cold nuclear matter as e ^ Ci. In the case of hot 
deconfined matter the light-cone loss per unit path is obtained by generalizing the expression in cold 
matter as 

Anas f dy-{n\e-^^[id-A''+{y-)A''^{0)]\n) 
euot = ^ . (87) 



' c 



Where \n) is a state in the thermal ensemble. In writing the above equation, we have made use of the 
boost invariance of the derived expressions along the direction of the jet. Unlike the expression for Ci 
which is derived in the Breit frame, Eq. (187|) is intended in the rest frame of the medium. In this frame. 
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Figure 11: The drag experienced by a hard quark jet in an HTL plasma. Figure adopted from Ref. [88] . 

one relates the light cone loss to the energy lost by a hard jet per unit length, for a jet traveling in the 
—z direction 



dq- dE^\dj)^\ dE{\^ 



dL~ dt + \dz\ \dz\{\v\~ 



dE 



+ 1) d\z\ 



(88) 



The other transport coefficient of importance to the light-cone or longitudinal propagation of a hard 
jet without radiation is obtained from the double derivative of the (5-function in Eq. (l85l) . This yields 
the diffusion in the light cone momentum per unit light cone path, denoted as the coefficient 62: 



62 



Ana J dy- {n\e-^^[F^~+ F^^+]\n) 



2Nr 



^9) 



Calculating the change in the light-cone fraction of the leading parton due to only these two coefficients 
amounts to incorporating a drag and a diffusion term in the energy distribution of the hard parton. 
This is often referred to as the Gaussian approximation with the location of the peak of the Gaussian 
given by e and the width given by 62- 

While operator products, in general, are evaluated at the hard scale at which they interact with the 
jet, they may be calculated at the scale of the medium if there exists a well defined microscopic model 
of the medium. These may then be evolved up to the hard scale of the jet. For example, for an HTL 
plasma, the elastic loss coefficient may be evaluated as 



d\k\^ f dx\kW-x){l-x^){N^ -1) 
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2n 



4iV. 



pAlklk'^) - PL{\klk')\ {1 +nB{k')} 



(90) 

-kx 



The equation above is obtained by evaluating Eq. (l87j) in the HTL limit. The factors pi and pt are the 
longitudinal and transverse HTL spectral densities (see Ref. [88] for details). The limits of the integral 
Qmax = V ET is set identically as for the calculation of q in an HTL medium. The elastic energy 
loss experienced by a hard jet with energy E propagating through an HTL plasma at a variety of 
temperatures and couplings are presented in Fig. [TTl 

With the discussion of elastic energy loss we complete the theoretical section of this review. In the 
remaining, we will outline the phenomenological modeling of the media encountered in jet modification 
experiments and review the main experimental measurements. 
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5 Parton energy loss in collisions: modeling of the medium 



Experimental measurements of parton energy loss are performed in cold nuclear matter [eA and pA 
collisions) and in hot nuclear matter produced in heavy ion collisions. In the experiment, the matter 
density is not uniform. In heavy ion collisions, there is a strong longitudinal expansion which causes 
a rapid decrease of the medium density. In addition, the point of origin of the hard parton within 
the matter is not known and therefore has to be integrated over. Due to these effects, it is necessary 
to include some form of averaging over the collision geometry in parton energy loss calculations when 
comparing to experimental results. 

We will first discuss the relevant aspects of the collision geometry, before turning to a more detailed 
discussion of the geometry averages used in the various calculations. 



5.1 Nuclear overlap geometry 

In central high-energy nuclear collisions, all of the incoming nucleons traverse the entire nucleus that 
moves in the opposite direction. The collision can be conveniently modeled in the approximation where 
the incoming nucleons travel along straight lines. 

The density distribution of nucleons in the nucleus follows the Woods-Saxon density profile 

p(r) = 



exp no- 



where R is the radius of the nucleus and d is the 'skin depth'; po is a normalization constant, such 
that the integrated density is equal to the atomic number A. The density in the transverse plane, the 
thickness function T, is then 



T(s) = J dzp{z, 



where z is along the longitudinal direction and s is a vector in the transverse plane. 

In a collision of two nuclei, the transverse distance between the centers of the nuclei is called the 
impact parameter b. Using the transverse impact vector b, we define the thickness function for the 
overlapping nuclei Tab and the collision density pcoii- 

T^b(s) = r^(s - ib) xTb{s+ ib) and p,M = Tab a. (91) 



Another measure of the overlap is the participant density 

Ppartis) = Ta{s - ib) • (^1 - e-^-(^+i^) '^^ + Tb{s + ^b) ■ (^1 - e"^-^^-^) '^^ . (92) 

Ta and Tg are thickness functions of the two colliding nuclei and a is the cross section for the process 
one is interested in. By integrating these two densities over the transverse coordinate s and taking a 
as the total inelastic cross section we obtain the total number of binary (inelastic) collisions A'^coii(^) 
and the total number of participants A'part(^), the number of nucleons that had at least one (inelastic) 
collision at a given impact parameter b. 

The number of participants and the number of binary collisions are the limiting cases of scaling of 
the particle production in heavy ion collisions, in a way that is related to the production cross section. 

• Participant scaling occurs for large-cross section processes, where every nucleon can produce par- 
ticles at most once. This is also called wounded nucleon scaling [90]. The multiplicity distribution 
in Au+Au collisions at RHIC follows this limit rather closely. 
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• Collision scaling occurs for small cross section processes, where every individual encounter between 
two nucleons contributes the same probability to the total production cross section. High-p^^ 
production of non-interacting probes at RHIC (photons) follows this limit (see Section [6.4. ip . 

In parton energy loss models, Ppart or pcou can be used as a first approximation for the medium 
density profile. Since the soft particle production scales more closely with the number of participants, 
Ppart is the most intuitive choice for a medium density profile. The distribution of initial hard scatterings, 
which are the points of origin for the partons that traverse the medium, is expected to follow pcou- 

In the above, we have used a continuum approach, the 'optical limit', where the nuclei are modeled as 
smooth density profiles. Alternatively, one can model the collisions with localized nucleons leading to a 
non-uniform density. This leads to important fluctuations of the density profile in peripheral collisions. 
For more details we refer the reader to [91] and references therein. 

The density profiles discussed above are suitable descriptions of the initial state geometry. The 
initial state density profile has pressure gradients that will lead to expansion, both in the transverse 
and longitudinal direction. In an expanding medium, the density decreases with time. The dominant 
effect is longitudinal expansion, which results in a density decreasing approximately as 1/r, where r is 
the proper time. Transverse expansion also leads to a reduction of the local density, but at same time 
increases the system size, so that the overall effect on energy loss observables is smaller. Transverse and 
longitudinal expansion can be modeled using hydrodynamic evolution. 



5.2 Modeling the produced matter: hydrodynamic evolution 



Relativistic fluid dynamics at low viscosity has emerged as the leading theoretical setup with which 
to describe the space-time evolution of the soft produced matter at RHIC. Global fits of a number of 
simulations to the experimental data have placed an upper bound on ratio of the viscosity t] to the 
entropy density s to be at most five times the absolute minimum allowed from quantum mechanics and 
predicted by the AdS/CFT conjecture [5^ : 

1 ?7 5 , , 

A-K s A-K 

Even ideal hydrodynamics where ?7 = yields a more or less consistent description of the soft spec- 
trum [93]. In these simulations, one solves the local differential equation, 

9^T'^'^(x,2/,z,t) = 0, (94) 

where, T^'^ is the energy momentum tensor, for the ideal case: 

T^"^ = {e + p) u^^u" - pg'^r (95) 

In the equation above, e is the local energy density, p is the local pressure and u is the fluid four velocity. 
The energy density is related to the pressure by the equation of state which is usually obtained by 
parametrizing lattice results. The set of equations are solved with the additional constraint of local 
baryon number conservation, 

d^^iuBU^') = 0. (96) 

A solution to the equations also requires a starting proper time Tq. Fits to data have placed this 
around tq = 0.6 fm/c. One also requires an initial configuration for the energy density and the baryon 
number density as a function of both the transverse location as well as the longitudinal location or 
rapidity. In the calculations presented in this review, e and ub were assumed to scale with each other 
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with an overall normalization fit from experimental data. For the energy density, one assumes the 
impact parameter dependent profile. 

e(^o, s, 7]) = eMAxp{s; b)H{ri). (97) 

Based on the equation of state and lattice data, the initial energy density can be related to the entropy 
density which may then be inferred from the produced multiplicity. The normalization of the the Baryon 
number density is also obtained by fitting to the final experimental data on p/n ratios. The transverse 
profile W{x, y; b) is obtained as a combination of the wounded nucleon density and the binary collision 
density. In 2+1 D hydrodynamics, the rapidity profile is assumed to be a constant while in 3+1 D 
simulations it is assumed to be a broad Gaussian, given by the form H{ri) = exp (— (?7)^/(2cr^)) with 
the width set from experimental data. 

The solution of the fluid dynamical equations results in a space-time profile of the energy momentum 
tensor and, as a result, of the energy density and pressure. As the temperature drops below the transition 
temperature, one enters the hadronic phase, at this point either the fluid dynamical simulation is 
continued with the hadron gas equation of state or one may switch to a hadronic cascade such as 
URQMD [9l]. Eventually as the temperature drops below the kinetic freezeout temperature the 
hadrons have to be decoupled from the calculation. In the case where hydrodynamics is continued into 
the hadronic phase, one uses a Cooper-Frye prescription [95] to convert the energy and momentum 
density in a unit cell into a distribution of free hadrons. 

The various input parameters along with the initial thermalization time Tq and final freezeout time 
are dialed to obtain the best fit with experimental data on soft hadrons of various flavors and from 
various centralities. Once all these parameters have been set, the space-time profile of the energy density 
can be used to calculate the quenching of hard jets. 



5.3 Energy loss in a non- uniform medium 



Hard partons are produced at all locations in the transverse plane, with a probability distribution 
following pcoii and no preference in azimuthal angle, and then propagate through the dense matter 
undergoing energy loss. In order to take into account the non-uniform, dynamic medium, one has to 
integrate over all production points and azimuthal angles, taking into account the medium density 
profile that the partons experience along their path out of the medium. 

All energy loss formalisms contain an integral over the scattering centers that the partons encounter 
along their paths. To obtain analytical results, the medium density profile is assumed to be either 
uniform or of a specific analytical shape {e.g. exponential decay in GLV). For a non-uniform medium, 
one can evaluate the full integral over the non-uniform medium density for every parton. This is what 
was done in the case of the HT formalism in [96] . 

For the ASW and GLV approaches, analytical results have been derived for a homogeneous expanding 
medium, where the density falls as 1/r along the path length [301 Sll EZ] • It has been shown in that in 
the infinite-energy regime those analytical results are well approximated by 

AE (X J dxxq{x), (98) 

where q{x) is a local transport coefficient and x is a coordinate along the path of the parton. Note that 
the integral is equal to ^qL"^ for a uniform medium. In addition to the scaling of the average energy 
loss, it has also been shown that in the homogeneous expanding medium, the differential gluon emission 
spectrum can be approximated using the calculation for a homogeneous static medium with effective 
transport coefficient 
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XQ+L 



gcff = 7 / dxxq{x). 

^ Jxo 



46 



This is often referred to as the 'dynamical scahng law'. A more general version of this scaling law was 
recently derived in ^5] . 

To calculate energy loss in a non-uniform, dynamical medium using the ASW formalism, it is 
assumed that the dynamical scaling law can be extended to this scenario. The first example of such 
a calculation is the Parton Quenching Model (PQM) by Dainese, Loizides and Paic [99], where pcou 
was used based on saturation model arguments |100j . More recent calculations use a hydrodynamical 
medium and where the local q{x) is proportional to the temperature T^, energy density e^^^ or the 
entropy density s [96l [Ml [102]. 

The AMY energy loss calculation is carried out in the finite temperature field theory approach 
and therefore directly relates the energy loss to the temperature T. In a hydrodynamical evolution, 
the local temperature T is known. The rates for a hard parton to decay into a quark or a gluon are 
then calculated using the local temperatures and the change in the distribution is computed using the 
Fokker-Planck equation over all paths [5S|, [HS] • 

5.4 Calculating observables 

Now we have all the ingredients for a calculation of parton energy loss that can be directly compared 
to experiment. For example, to calculate the nuclear modification factor Raa-, one starts from the 
factorized calculation in the vacuum, Eq. Parton energy loss is then included by substituting the 
vacuum fragmentation function with the medium-modified fragmentation function, which is given in 
Section [3l averaged over the medium density profile, as explained in the previous section. In general, 
the average over the medium density profile involves integrating out the parton production point and 
direction. 

In the next Chapter, we will compare energy loss calculations to experimental data on the nuclear 
modification factor Raa for light and heavy hadrons, and di-hadron and 7-hadron measurements. The 
di-hadron measurements involved pairs of partons, which are produced at a single point in the medium 
and then propagate outward in opposite directions in the transverse plane (back-to-back). As a result, 
the pathlength and density profile that both partons see are not independent. For example, if the pair 
originates from the periphery of the collisions zone and travels outward, one parton is likely to see a 
long pathlength in the dense medium, while the other will see a short pathlength. These correlations 
have to be included in the model calculations as well. 

Early comparisons of energy loss calculations to experimental data were based on a simplified 
medium geometry, using hard spheres or cylinders [103[ 11041 1105] , instead of the Woods-Saxon geometry 
in Eq. flS.ip and no hydrodynamical evolution. There are also calculations using a Woods-Saxon overlap 
geometry, without hydrodynamical evolution [211 [TS]. These are more realistic than the hard-sphere 
geometries in the sense that they include the low-density regions at the edge of the nuclei. In this 
review, we will focus on calculations which sample the space-time evolution using a hydrodynamical 
calculation for the medium density evolution [101[ [96l 1102] . 
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6 Comparison to experiment 



6.1 Choosing observables 

In the following, we review a selected set of experimental data and compare to energy loss calculations 
to show how the measurements constrain the energy loss formalisms and the properties of the medium. 

A wealth of different measurements have been performed to provide insight in the energy loss 
mechanism. In general, care should be taken to select observables that are sensitive to the process that 
one wants to investigate. This selection process requires a careful exchange of ideas between theory and 
experiment. For example, the first measurements at RHIC of the nuclear modification factor [21 Ej and 
di-hadron correlations (disappearance of the away-side jet) |106] have clearly shown that energy loss is 
large. It was subsequently realized, however, that a variety of models can describe the data, because 
these measurements are not very sensitive to the details of the energy loss distribution. It has been 
shown that it is in practice impossible to extract the energy loss probability distribution from such 
measurements [1U7] even if the experimental uncertainties would be exceedingly small. It is therefore 
key to identify observables that are sensitive to specific aspects of the energy loss process. For example, 
it has been argued that comparing the away-side suppression to the single hadron suppression allows 
to determine the path-length dependence of energy loss (L vs L^) |108] . 

In Figure [10] it was shown that the different formalisms for parton energy loss calculations predict 
different medium modified fragmentation functions. Such differences can only be determined from 
experiments if the initial parton energy can be controlled or measured. This technique has been used 
successfully in Deeply Inelastic Scattering experiments off nuclei to measure parton energy loss in cold 
nuclear matter (see Section [6.3. ip . In heavy ion collisions, there are two basic ways in which one can 
control the initial parton energy. One is to use 7-jet events (or Z°-jet at LHC) where the transverse 
momentum of the photon and the initial parton (jet) are equal (as long as initial state broadening 
and higher order large angle radiations are negligible). The alternative is to reconstruct jets in the 
final state. To what extent the reconstructed jet energy is equal to the initial parton energy depends 
both on the jet reconstruction algorithm and on the effect of medium modifications on the jet. Several 
measurements are likely needed to disentangle various effects (jet spectra, jet fragmentation functions 
and away-side jet rates have been discussed). 

The primary consideration when selecting observables should be the intrinsic sensitivity: the larger 
the effect of a given model aspect on the observed quantity, the more promising the observable is. This 
type of optimization is performed based on theoretical insight and model calculations. 

A secondary consideration is the experimental precision. Experimental uncertainties vary from ob- 
servable to observable, depending mainly on the cross section and the amount of background. Experience 
at RHIC shows that for the least complicated measurement, such as pion spectra, an experimental un- 
certainty in the 5-10% range can be reached. More complicated measurements, with larger backgrounds 
that cannot always be accurately modeled in heavy ion events, have correspondingly larger uncertainties. 
Further reduction of the uncertainties often requires a large effort, often requiring dedicated detector 
construction and calibration for a specific measurement, and is therefore only pursued by experiments 
if there are compelling theoretical reasons to do so {e.g. precision tests of the electroweak sector of the 
standard model). A few specific examples of how the experimental precision may evolve at RHIC are 
given in Section 16.51 

6.2 Perturbative QCD in p+p collisions 

Factorized perturbative QCD calculations have been compared extensively to experimental data on p-p 
collisions at energies ranging from ISR ^/s ~ 60 GeV to Tevatron energies ^/s = 1.9 TeV |113[ 1114] I115[ 
I116j . RHIC measurements at mid-rapidity generally agree well with perturbative QCD calculations. 
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Figure 12: Left panel: Comparison of the 7r° pT-spectrum in p+p collisions at y/s = 200 GeV as 
measured by PHENIX [109] and charged pions measured by STAR [llOj to NLO pQCD calculations 
[TT] at three scales = pt (solid line) and /i = 0.5pT and /i = 2pT (dashed lines). The uncertainty bars 
on the left side of the lower panel show the overall normalization uncertainty in the STAR (dark grey 
band) and PHENIX (light grey band) measurements. Right panel: Comparison of the pr-spectrum 
of non-photonic electrons as measured PHENIX to a NLO pQCD calculation with resummed 

Next-to-leading logarithms (FONLL) |112] . The dashed lines indicate the uncertainty on the FONLL 
calculation, which is estimated by varying the factorization scale and the quark masses. The lower 
panels show the relative difference with the central curves. The normalization uncertainty is indicated 
in the lower panel on the right. 
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Two examples are shown in this section. 

Figure [12] (left panel) shows a comparison of pion spectra measured by STAR and PHENIX, with 
STAR measuring charged pions in a Time Projection Chamber (TPC) |llUj and PHENIX measuring 
neutral pions with an Electromagnetic Calorimeter (EMCal) |109] . The measurements are compared 
to Next-to-Leading Order (NLO) pQCD calculations by Jager et al. [11], using the CTEQ6M parton 
density functions [9] and fragmentation functions by Kniehl, Kramer and Potter (KKP) |117j . The 
measurements by STAR and PHENIX agree well with each other and with the NLO pQCD calculation. 
The uncertainties on the calculation are estimated by varying the factorization and renormalization 
scales between /i = O.Spr and /i = 2pT and are indicated by dashed lines on the figure. 

Heavy quark production is measured at RHIC via the semi-leptonic decays D ^ e and B ^ e. The 
inclusive electron distributions that are measured in the experiment contain substantial contributions 
from photon conversions and Dalitz decays of 7r° and t] mesons which are subtracted to obtain the 'non- 
photonic' electron yield. The p^-dependent inclusive cross section of non-photonic electron production 
as measured by PHENIX |lllj in p+p collisions at a/s = 200 GeV is shown in the right panel of Fig. [12] 
and compared to NLO pQCD calculations including resummed Next-to-Leading Logarithms (FONLL) 
[112] . A similar measurement was performed by the STAR experiment |118] and showed a larger yield 
of non- photonic electrons. This measurement, however, has recently been superseded by a new, still 
preliminary, result |119] . The uncertainties on the theory expectation are estimated by varying the 
renormalization and factorization scale and the quark masses. The yields measured by PHENIX are 
larger than the nominal theory expectation, but close to the upper bound on the expectation. 

6.3 Cold nuclear matter effects 

Before turning to nucleus-nucleus collisions to study parton energy loss in hot nuclear matter, we will 
briefly discuss the relevant effects in cold nuclear matter. We distinguish initial state and final state 
effects in cold nuclear matter. The dominant initial state effect is the modification of the parton density 
distribution in the nucleons, or the 'EMC-effect', named after the first experiment that measured this 
effect |120j . The most important final state effect is parton energy loss in cold nuclear matter, which 
has been found to be significant in the HERMES experiment [31]. Both effects are discussed below. 

For completeness, we also mention here the 'Cronin effect', a possibly mass- dependent enhancement 
of hadron production at intermediate pt in proton-nucleus collisions |121] which is generally attributed 
to transverse momentum broadening due to multiple scattering of partons on the nuclear target. For a 
recent review, see |122] . 

6.3.1 Energy loss in cold nuclear matter 

Figure [13] shows a comparison the measured suppression of pions from fragmentation in cold nuclear 
matter in eA collisions by the HERMES experiment [31] with calculations using the BDMPS-ASW 
quenching weights (left panel) and the higher twist formalism (right panel). In deeply inelastic scattering 
experiments, the jet energy can be determined from the momentum vector of the scattered electron; the 
measurement is presented in terms of the momentum fraction z = Eh/v where is the hadron energy 
and V the energy of the virtual photon. The extracted transport coefficient for the BDMPS quenching 
weights is g = 0.6 GeV^/fm in the center of the nucleus. Using the higher twist formalism, including the 
effect of multiple emissions, a quark q ~ 0.08 GeV^/fm is obtained. In an earlier calculation, involving 
only one emission in medium |103] . a g = 0.12 GeV^/fm was obtained |124] . 

6.3.2 Experimental study of initial state: d-|-Au 

At RHIC, d+Au collisions are used to measure the effect of possible initial state effects on the pion and 
photon spectra. The left panel of Fig. [14] shows recent measurements of the nuclear modification factor 
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Figure 13: Comparison of the measured suppression of pion (full symbols vr"'", open symbols tt") pro- 
duction in eA collisions by HERMES [HT] with calculations using ASW-BDMPS quenching weights (left 
panel) and the higher twist formalism (right panel). Figures from [TO] (©Societa Italiana di Fisica) 
and [T23] . 
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Figure 14: Nuclear modification factor RdAu for pions (left panel) and photons (right panel) in d+Au 
collision at ^/sNN = 200 GeV, compared to calculations with nuclear parton density functions (left 
panel) and various other cold nuclear matter effects (right panel). Figures reproduced with permission 
from [125] and |126] (right panel ©American Physical Society). 
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RdAu for 7r° from PHENIX |127] and charged pions from STAR |128j and compares to calculations using 
a recent analysis of nuclear parton density functions (nPDF) (EPS09 |125j and HKN07 |129] ). It should 
be noted that the DIS results do not fully constrain the nPDFs and the RdAu measurement has been 
included in the determination of the nPDFs. The good agreement with the data therefore indicates that 
the measured values of RdAu are consistent with the nuclear modification of PDFs which is measured 
in DIS off nuclei |120] . No significant final state interaction is apparent in the data. 

The apparent absence of cold nuclear matter energy loss in the RdAu in Fig- [H at first sight is 
surprising, given the large suppression of hadron production at high momentum fraction z in the DIS 
data in Fig. [131 However, the d + Au results are at ?/ = in the laboratory; in order to compare with 
the HERMES results, we boost to the rest frame of the nucleus where the jet with only a transverse 
momentum pj- at midrapidity has an extremely large energy due to the boost. The HERMES measure- 
ments show that the effect of energy loss reduces with jet energy [31] and is expected to vanish at the 
very high energy of the boosted jet. A detailed modeling of this effect is given in |122] . 

The right panel of Fig. [H] shows a comparison of the measured nuclear modification factor RdAu 
for photons in d+Au collisions at RHIC {^/sNN = 200 GeV) with calculations including various cold 
nuclear matter effects |126] . The calculations show a modest enhancement aX pt < 5 GeV/c due to 
the Cronin effect and modest suppression at higher px due to a combination of isospin effects and 
shadowing. The curves with energy loss (labeled AE{IS)) refer to energy loss of the incoming parton 
as it propagates through the nuclear matter before the hard scattering |130] . The total effect on the 
photon spectrum is around 20%, which is slightly smaller than the current experimental uncertainties 
in the measurement. 

6.4 Heavy ion collisions at RHIC 

At RHIC, several types of measurements that are sensitive to parton energy loss in hot and dense QCD 
matter have been performed. In the following, we will discuss the results on inclusive hadron suppression, 
di-hadron suppression as well as results for heavy quarks. Comparisons to model calculations of energy 
loss are made where possible, with emphasis on calculations based on a full hydrodynamic evolution of 
the density and radiative energy loss according to the four formalisms discussed in Chapter |2l 

We also focus on the region of high px, where perturbative QCD is most likely to be applicable. 
Experimental results at RHIC show anomalously large baryon to meson ratios at ^ 6 GeV/c |13H 
I132[ 1133] , indicating that there are important contributions to hadron production beyond factorized 
perturbative QCD in that regime. We therefore limit the discussion to Pt > Q GeV/c where possible. 

6.4.1 Verifying initial production rates: Photon production 

Figure [15] shows the centrality dependence of the nuclear modification factor Raa (Eq- [TT]) for prompt 
photons in central Au+Au collisions as measured by PHENIX at ^/sNN = 200 GeV, compared to 
various model curves. The measured Raa for photons is close to unity, as expected for hard probes 
{Ncoii scaling) without energy loss. 

The model curves in Fig. [15] show the effect of initial state energy loss (left panel) and final state 
energy loss for fragmentation photons (right panel) and medium-induced photons, including medium 
induced QED bremsstrahlung and parton-photon conversion in scattering (right panel). These mech- 
anisms, together with isospin effects, are expected to lead to a suppression of the photon yields by up 
to 25% at high px > 10 GeV, with the exact value depending on details of the modeling. The model 
predictions are in agreement with the experimental data, but the current large uncertainties on the 
measurement preclude quantitative conclusions about the contributions of the various effects. 
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Figure 15: Nuclear modification factor Raa for prompt photons in central Au+Au collisions at ^/snn = 
200 GeV as a function oi pt, |134] compared to various model calculations. Left panel: Model curves by 
Vitev, including initial state energy loss |126] . Right panel: Model curves by Qin et al showing separate 
contributions of direct, fragmentation and medium induced photons [30]. The black solid line shows 
the effect of adjusting the p + p reference in the calculation to the measured spectrum. Both panels 
©American Physical Society. 
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Table 1: Values of the initial transport coefficient go {t = tq = 0.6 fm/c) at the center of the most central 
collision (0-5%) as determined by comparing calculations of Raa with a medium density space-time 
profile taken from a 3D hydrodynamical evolution and three different parton energy loss models [96] . 



6.4.2 Light hadron suppression 

Parton energy loss leads to a suppression of the pion yield at high pt- A recent measurement of Raa 
of 71^ by the PHENIX experiment in shown in Fig. [161 A clear suppression is visible at high pt- The 
suppression increases with centrality and is approximately pT-independent for > 4 GeV/c. The value 
of Raa at high pt for central collisions is around 0.2. 

Figure [17] compares calculations of the nuclear modification factor using a 3D hydrodynamical 
evolution for the medium density and three different parton energy loss formalisms to experimental 
data on 7r° suppression from PHENIX in Au-|-Au collisions at ^/sN^ = 200 GeV and two different 
centralities. The hydrodynamical simulation is tuned to fit RHIC data with a < 2 GeV and 
provides space-time profiles for the medium properties. The only free parameter in the calculation is an 
overall proportionality constant between the local transport coefficient and the medium density. This 
parameter was set using the value of Raa at pr = 9.5 GeV in central collisions. The behaviour with pt 
and centrality are parameter-free predictions. It is interesting to note that all three model calculations 
show only a very weak dependence of Raa on px, in agreement with the data. In the model calculations, 
the independence of Raa on px does not have a simple interpretation; it comes about as a result of the 
interplay between the shape of the parton spectrum and the energy-dependence of energy loss, which 
is partially due to the kinematic cut-off AE < E. 

The resulting values of the transport coefficient at thermalization (r = 0.6 fm/c) at the center 
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Figure 16: Nuclear modification factor Raa as a function of for ir^ as measured by PHENIX in 
Au+Au collisions at ^/sNN = 200 GeV. The different panels show different centrality selections. The 
error bands around the dashed line at Raa = 1 indicate the uncertainty from the number of binary 
collisions (left band) and the normalization uncertainty in the p + p reference data (right band). Figure 
reproduced with permission from |135] (©American Physical Society). 
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Figure 17: Nuclear modification factor Raa as a function of as measured by PHENIX for tt compared 
to calculations using a 3D hydrodynamical medium evolution and 3 different energy loss formalisms: 
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of the most central collisions (0-5% centrality) are given in Table [TJ The extracted values for q depend 
on whether the local transport coefficient is assumed to be proportional to the energy density e (to the 
appropriate power 3/4), the third power of the temperature T or the entropy density s. In addition, 
there are significant differences between the energy loss formalisms: the Hard Thermal Loop formalism 
(AMY) gives go = 4.1 GeV^/fm, while the higher twist formalism gives a lower value of 2.3 GeV^/fm 
and the multiple soft scattering approach (ASW) gives a much larger value of 10 GeV^/fm. 

It is interesting to note that in this comparison, the HT and AMY formalisms give similar medium 
density parameters, while it can be seen in Fig. [10] that AMY generally has a larger suppression at given 
T and L than HT. Also, in Fig. [TOl the suppression for ASW is closer to HT for L = 2 than to AMY, 
while in this comparison ASW needs a much larger value of go- The qualitative difference between 
the comparison to Raa in Fig- CZl and the fragmentation functions in the Fig. [10] is probably due to 
differences in the treatment of the medium in the hydrodjTiamical calculations of Fig. [T7] and Ref . [96] : 
for the AMY calculation, no energy loss was taken into account for r < tq and in the hadronic phase 
at late times, while for the HT calculation in the same paper, the density was taken to be constant for 
r < To and energy loss was also calculated in the hadronic phase. The ASW calculation in that paper 
calculated no energy loss for t < tq but did continue the energy loss calculation in the hadronic phase. 

The treatment of the density for early times r < tq, when the medium is very dense, has a large 
impact on the extracted medium density. A factor 2 difference was found between the case where g = 
is taken for r < tq compared to constant g = go for the initial stage in Ref. |102] . The treatment of 
late times may also have a significant quantitative impact; if energy loss continues after the medium 
has hadronized. Due to the rarefaction wave in hydrodynamic simulations, the path length does not 
increase much beyond the point of hadronization, however contributions from the hadronic phase may 
still have a noticeable effect on the energy loss. A somewhat simple-minded study of this effect was 
carried out in Ref. |136] . where a rescaling parameter called chg was used to dial the effective g in the 
hadronic phase compared to the g in the partonic phase at the same temperature. It was found that 
choosing chg to be vanishing required almost a doubling of the effective g in the partonic phase for the 
results to remain consistent with the data. An alternative approach has been adopted by the authors 
of Ref. |137] . who used the g extracted from jet modification in DIS on a large nucleus to estimate the 
g in the hadronic phase of a heavy-ion collision. 

Thus far, the importance of the early and late time dynamics in the modeling of parton energy loss 
has not been extensively studied in the literature. Future studies will need to explore this aspect further, 
for example by using the exact same treatment for side-by-side comparisons of other model aspects. 
Uncertainties in the initial density of the system should be taken into account when determining the 
medium density from energy loss measurements and will likely give a sizable contribution to the overall 
uncertainty. The present authors feel that taking g = for r < To is unrealistic, since it also implies 
a very sharp rise of the density at thermalization r ~ To. A more realistic lower limit of the medium 
density at early times is needed, but there currently exists no commonly accepted estimate. 

At the moment of writing, a comparable calculation of Raa with a hydrodynamic medium evolution 
using the GLV energy loss formalism was not available. Early work by Hirano and Nara |138] uses 
full hydrodynamic evolution, but only a schematic form of the GLV energy loss formalism, and can 
therefore not be compared directly. 

6.4.3 Di-hadron correlations and Iaa 

The measurement of inclusive spectra and the nuclear modification factor integrates over a large range 
of parton energies and over the entire collision geometry. Di-hadron correlation measurements provide 
more differential information on the partonic kinematics and a different geometrical bias. Comparing 
results may therefore provide more insight in for example the path length dependence of energy loss. 
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Figure 18: (Left) Distributions of azimuthal angle difference A0 between charged associated particles 
and a high-pT- charged trigger particle {8 < px < 15 GeV/c) for different system sizes (d+Au, mid- 
peripheral and central Au+Au) at y/s^^ = 200 GeV/c and different p^-ranges for the associated 
particles, including the uncorrelated background. The distributions are normalized per trigger particle. 
(Right) Associated yield on the near side yield |A0| < 0.63 (left panels) and recoil \A(j) — 7r| < 0.63 
(right panels), as a function of zt = for trigger particles with {8 < pr < 15 GeV/c) and three 

different system sizes at ^/JJvw = 200 GeV/c. The lower panels show the ratio of the Au+Au result to 
the d+Au reference measurement. Both panels are reprinted with permission from |139] (©American 
Physical Society). 



The left panel of Fig. [18] shows example distributions of the azimuthal angle difference A0 between 
charged assoczated particles and high pr {8 < pr < 15 GeV/c) trzf/^'er particles. In all panels, two peaks 
are visible, which is a clear di-jet signature. For the lower p^ associated hadrons (upper panel), a clear 
increase is visible in the combinatorial background when comparing d+Au (left column) and central 
Au+Au collisions (right column). The near-side peak around A0 = above background has similar 
magnitude in d+Au and central Au+Au, while a clear suppression of the yield is visible in the recoil 
at A(j) = 71. It is also interesting to note that the width of the correlation peaks is similar in d+Au 
and Au+Au; no sign is seen of the transverse momentum kicks from the medium, neither at the near 
side (jet broadening) nor in the recoil (acoplanarity). Since both jet broadening and acoplanarity are at 
some level unavoidable consequences of the energy loss process, this leads to the conclusion that either 
the transverse kicks are small compared to the trigger and associated hadron momenta, so that the 
angular broadening is small, or that the selection of high-momentum particles biases towards partons 
that had little or no interaction with the medium. Note that only 20-30% of the recoil yield is visible 
in Au+Au, so a significant fraction of back-to-back di-hadrons that is present in the d+Au reference 
measurement does not pass the p^-selection in Au+Au collisions due to energy loss. 

A more quantitative representation of the di-hadron results is shown in the right panel of Fig. [18] 
which shows the per-trigger associated yield, after subtraction of the combinatorial background, on the 
near side (A0 ~ 0, left) and recoil (A0 ~ vr, right) with a trigger particle of 8 < < 15 GeV/c for 
d+Au collisions and Au+Au collisions at two different centralities. The lower panels show the ratio of 
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the Au+Au results to the d+Au reference. One the near side, no modification of the associated yield 
is seen, while the recoil yield is suppressed by a factor ~ 4. 

It is important to realize that the near-side and recoil yield in this measurement carry different 
information due to the presence of a trigger particle. The near side measures the yield associated with a 
trigger particle in the same jet. The absence of suppression on the near side does not necessarily imply 
that the jets emerge without energy loss; it is also possible that partons lose energy and subsequently 
fragment in the vacuum, while the lost energy is carried largely by hadrons at lower pt or outside the 
angular region of the measurement. The trigger particle selection would then select partons that have 
a momentum distribution after energy loss that is similar to the d+Au measurement. 

Theoretical calculations for the near-side associated yield are quite involved and require the in- 
troduction of a new non-perturbative object, the dihadron fragmentation function |14m 1141] which 
contains the correlation between the trigger and the associated hadron. The scale evolution of the di- 
hadron fragmentation function mixes with the single hadron fragmentation functions. The calculation 
of medium-modified di-hadron fragmentation has only been carried out in the HT formalism, we refer 
the reader to Ref. |142] for details. Figure [191 compares theoretical calculations for the case of one 
scattering and one emission in the medium with experimental data on the centrality and associated pt 
dependence of the associated yield. Note that since the transport coefficients have been fit to the Raai 
there are no free parameters in this calculation. The agreement is excellent for the case of the highest 
associated px and gradually moves away from the data as the range of px for the associated hadron is 
lowered. 
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Figure 19: The associated yield per trigger for a trigger with 8 GeV< pr < 15 GeV for three different 
ranges of associated px and for all centralities. The dashed line represents the results for the associated 
yield of charged hadrons where the initial condition for the dihadron fragmentation function is taken 
from JETSET. The dot-dashed line is the result for charged pions where the initial condition is taken 
from JETSET. The red lines are the results for the associated yield of charged hadrons where the 
initial condition has been fitted to the associated yield of charged hadrons in p-p collisions in STAR. 
Reproduced from |142j (©American Physical Society). 



The recoil yield, on the other hand, consists of fragments from the recoil parton (back-to-back 
with the triggered parton) which propagate independently. Energy loss of the recoil parton must lead 
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Figure 20: Model curves from a model with a hydrodynamical evolution of the medium and energy loss 
using the multiple soft scattering approximation (ASW quenching weights) compared to measurements 
of 7r° Raa (left panel) and charged particle Iaa- The broad green line indicates the best simultaneous 
fit to both experimental results Figure reproduced with permission from |102j (©lOP Publishing). 

to reduced high-p^ recoil particle yield per trigger. At lower associated pr, enhanced yield has been 
observed, suggesting that the radiated energy remains correlated with the trigger particle |143l 1144] . 

It is also interesting to note that the away-side suppression in the di-hadron measurement is inde- 
pendent of zt in the range shown in Fig [181 This suggests that the longitudinal momentum distribution 
of high-pT fragments is unmodified in Au-|-Au collisions, which could mean that the measured yield 
is mostly from partons that have little or no interaction with the medium, while, on the other hand, 
partons that did interact with the medium only produce associated yield at lower px- However, it should 
be kept in mind that the steeply falling parton spectrum, combined with approximately exponential 
fragmentation functions imply that the trigger pt selection provides only limited control over the ini- 
tial state parton kinematics. More complete modeling is needed to provide quantitative limits on the 
qualitative trends discussed above. First measurements with experimentally controlled parton energies 
via 7-hadron correlations and direct jet reconstruction are now becoming available, which can be more 
straightforwardly interpreted in terms of parton energies. These will be discussed in more detail in 
Sections l6X8l 

6.4.4 Model comparisons of Raa and Iaa 

Early comparisons of theory to data |1U3[I1U4]I1U5] have focused on the nuclear modification factor Raa- 
As shown in Fig. [T71 most energy loss formalisms can be tuned to describe Raa-, including the observed 
independence of Raa on Pt- To further validate energy loss calculations, a simultaneous comparison to 
observables with different geometrical and fragmentation biases is needed. 

Figure [20] shows a comparison of a model using a hydrodynamical evolution for the medium density 
profile and the ASW multiple soft scattering formalism for gluon radiation |102] . The model can describe 
the measurements of the nuclear modification factor Raa and the recoil suppression Iaa with a common 
scale factor K between the energy density in the hydrodynamical model and the transport coefficient 
q. The obtained values for K are in the range 3-5, indicating that the transport coefficient is larger 
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Figure 21: Nuclear modification factor Raa for emitted at angles close to the reaction plane and 
perpendicular to the reaction plane (75° < (j) < 90°) as measured by PHENIX in 20-30% central 
Au+Au collisions at ^/snn = 200 GeV |146j . The curves in the left panel represent calculations with a 
hydrodynamic medium, comparing radiative energy loss with three scenarios for elastic energy loss (see 
also Fig. |22|) |108] . The curves in the right panel represent calculations with a hydrodynamic medium 
at fixed impact parameter {b = 7.5 fm) in the same centrality range and using three different energy 
loss formalisms PB] . 

than the perturbative estimate of Baier (g ^ 2e^^^) |145] . These values are in approximate agreement 
with the values of go in Tabled] and also with an earlier fit by Renk and Eskola |101] . 

6.4.5 Path length dependence of energy loss 

Conceptually, a measurement of the path-length dependence of energy loss is very interesting, since it 
would directly show to what extent the LPM interference is present. As discussed in Section I2l2| LPM 
interference effects lead to energy loss proportional to in the infinite energy limit. Elastic energy 
loss, on the other hand, would be proportional to L. 

Experimentally, the path length of the partons through the medium can be varied by selecting 
different centralities and/or using different size nuclei. For this purpose, Cu+Cu collisions have been 
used at RHIC [147^ 1148] , but so far no detailed comparisons with model calculations with a realistic 
hydrodynamic medium have been performed. 

More differential control of the geometry is obtained by measuring Raa <^ 

function of the emission 

angle with respect to the reaction plane, or by measuring elliptic flow V2 at high p^- The most accurate 
measurement of this type has been performed by PHENIX. Fig. [21] shows the measured Raa for 7r° 
produced in 15 degree angular intervals around the reaction plane (short path length, square markers) 
and perpendicular to the reaction plane (long path length, circle markers) in mid-central (20-30% 
centrality) collisions are compared to two sets of model calculations using a hydrodynamical medium . 
In the left panel, model calculations for radiative energy loss (solid curves) and elastic energy loss are 
shown |108] . The three different calculations using elastic energy all have energy loss proportional to 
L with a gaussian spread of different widths. The elastic models are seen to have in general a smaller 
difference between Raa for hadrons emitted in the reaction plane and perpendicular to the reaction 
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Figure 22: Comparison of energy loss calculations with the nuclear modification factor Raa from Fig. 
[T71 (left panel) and the recoil yield measurement from Fig. [TS] (right panel, yellow points). One of the 
calculations uses radiative energy loss and the other three use different schematic implementations of 
elastic energy loss. The calculations use the same geometry from a hydrodynamic evolution and have 
been tuned to reproduce the inclusive hadron suppression measurement Raa- Figures reprinted from 
[108] with permission (©American Physical Society). 

plane than the radiative calculation, as expected from the different L dependence. The data are more in 
agreement with the radiative energy loss curves at the highest although the differences between the 
models are small. The out-of-plane measurement shows significant deviations from the model curves at 
Pt < Q GeV/c, possibly due to contributions from hydrodynamical flow of the soft matter. 

The right panel of Figure [2T] shows a comparison to three different calculations with radiative energy 
loss (the same calculations as Fig. [T71) . It can be seen in the figure that the AS W energy loss formalism 
leads to the largest difference between in-plane and out-of-plane suppression, which is closest to the 
experimental result. The statistical uncertainties in the measurement are still sizable and will be reduced 
in the future when larger data samples are analyzed. 

According to model calculations, the recoil suppression Iaa also has a significant sensitivity to the 
path length dependence of energy loss. The reason for this sensitivity is that the trigger particle selects 
partons on the near side that have little or no energy loss and therefore a short path length. As a 
result, the recoil parton has a longer-than-average path length. This so-called 'surface bias' effect was 
first pointed out by Dainese, Loizides and Paic in [99]. Later, it was pointed out by Renk |149j that 
there is an important interplay between the expansion of the medium and this path- length bias effect. 
The longest path-lengths arise for partons that originate close to the surface, but then travel inwards. 
These partons only pass through the center of the collision after a few fm/c. In the presence of rapid 
longitudinal expansion, the medium will have cooled down significantly by that time and thus despite 
the long path length, the average medium density along the path is significantly reduced. 

Fig. [22] (right panel) shows a comparison of the measured recoil yield with high px trigger particles 
8 < Pt < 15 GeV/c with calculations using ASW quenching weights for radiative energy loss and the 
same three implementations of elastic energy loss that were used in the left panel of Fig. [211 based on 
Gaussian energy loss fiuctuations with different mean and width for the Gaussian. In each calculation, 
the medium density has been tuned to reproduce the inclusive hadron suppression Raa (left panel). 
The figure (right panel) clearly shows that all three elastic energy loss models under predict the recoil 
suppression, due to the combination of the path length bias and the linear dependence of elastic energy 
loss on L, thus indicating that radiative energy loss, with a quadratic dependence on L, is dominant 
for light hadrons. Comparing the left panel of Fig. [21] and the right panel of Fig. [22] we also note that 
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the recoil suppression Iaa is more sensitive to the path length dependence of energy loss than the angle 
dependence of Raa- 

6.4.6 Colour factor dependence: quark versus gluon energy loss 

Radiative energy loss is expected to be Ca/Cf = 9/4 times larger for high energy gluons than high- 
energy light quarks. I has been suggested that this difference is experimentally accessible by comparing 
nuclear modification factors for (anti-)protons and pions. In e^e~ collisions at LEP, it has been observed 
that gluon jets are more likely to produce (anti-)protons than quark jets |15U] . One could therefore use 
anti-protons to 'tag' gluon jets and expects to see a larger suppression for (anti-)protons than for pions. 

There are several caveats to this simple expectation. First of all, the argument crucially depends 
on the difference between (anti-)proton production from gluon and quark fragmentation. The gluon 
to proton fragmentation function is not accurately know; in fact, recent extractions of fragmentation 
functions from experimental data use the RHIC p-p data to constrain this [15H 1152] . It has also been 
shown in more detailed calculations using the AKK fragmentation functions, that while the expected 
difference between quarks and gluon Raa is almost a factor 2, the effect is diluted when comparing 
protons and pions, because pions are produced by gluon as well as quark fragmentation |153] . Finally, 
it should be realised that quarks showers contain gluons and vice versa, so that the distinction between 
quark and gluon energy loss in a full shower evolution not be clear. One specific implementation of this 
effect are 'jet conversion' processes in the medium. A calculation of this effect shows a reduced p/vr 
ratio in heavy ion collisions compared to p-p collisions |154] . Measurements performed by STAR |132] , 
do not show a larger suppression of (anti-)protons compared to pions. 

6.4.7 Heavy quark energy loss 

Due to their large mass (m >> Aqcd)? heavy quarks move at speeds smaller than the speed of 
light through the medium. As a result, they cannot to emit gluons at small angles, the so-called 'dead 
cone effect', and therefore lose less energy than light quarks (see Section WA\ |85] . To verify whether 
radiative energy loss is dominant for heavy quarks, measurements have been performed of heavy flavor 
production using semi-leptonic decays into electrons. 

The upper panel of Figure |2S] shows a comparison of predictions using only radiative energy loss by 
Armesto et al (solid line in Fig. [2^ . using the parameters obtained from the light quarks measurements 
in Fig. [201 to measurements of the nuclear modification factor of non-photonic electrons (electrons 
which are not from photon conversion or light hadron decay) from heavy flavor decays by PHENIX 
at RHIC |102] . The measured suppression of non-photonic electrons is larger than the model predicts. 
The extracted i^'-factor would be larger than 25, while the value for light quarks was between 2 and 5. 
These results indicate that the suppression of heavy quarks is larger than expected, possibly similar to 
light quark suppression. It is worth noting that the systematic uncertainties on the measurements are 
sizable. These uncertainties are expected to be quite strongly correlated point-to-point. As a result, the 
data are not strictly inconsistent with the expectation from radiative energy loss. For a more detailed 
discussion of the systematic uncertainties and the implications for energy loss we refer to |102] . 

Various mechanisms have been proposed to explain the observed large suppression. For example, 
it has been argued that elastic energy loss may play a significant role for heavy quarks, because the 
radiative energy loss is reduced [7S]. A detailed evaluation of the combined effect of radiative and elastic 
energy loss has been performed using the opacity expansion for radiative energy loss and elastic energy 
loss based on the HTL calculation by Thoma, Gyulassy (TG) [82] and Braaten, Thoma (BT) [83]. The 
same mechanism should be present for light quarks, so the model was fitted first to the light quarks 
and the predicted suppression of non-photonic electrons is compared to the measurements in the top 
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Figure 23: Nuclear modification factor Raa for non-photonic electrons from heavy flavor decays as 
measured by PHENIX |155] in central Au+Au collisions at = 200 GeV compared to three 

different model calculations. The uncertainty bands on the lower left indicate the overall normalisation 
uncertainty. Top panel: radiative energy loss only (see also Fig. I2U]1 |102] and radiative and elastic 
energy loss, using two different elastic energy loss models [TS]. Bottom panel: radiative and elastic 
energy loss [87], using transport coefficients from HTL effective theory (solid line) and with an alterative 
prescription (dashed line, see text). 
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panel of Fig. |23] (dashed and dotted lines). The discrepancy between the model and the experiment 
is smaller than when only radiative energy loss is taken into account (compare to solid line), but the 
calculation still gives less suppression than the measurement. 

In the inclusion of elastic energy loss in the calculation above one assumes not only that the jet is 
weakly coupled to the medium, but also that the medium may be described as a leading order HTL 
plasma. The expressions for the elastic energy loss in both the BT and TG calculations in the upper 
panel of Fig. [23] assume that the medium may be described as weakly interacting quasiparticles which 
interact directly with the jet leading to the elastic energy loss. In this limit, all the transport coefficients 
have a characteristic dependence on the logarithm of the parton energy, see Eq. f l73|) . which reduces the 
energy loss for lower jet energies. Since all the relevant transport coefficients have the same logarithm, 
there remains only one fit parameter which is the value of as in the medium. This is set by fitting to 
one data point in the light flavor sector. 

The ascribed limits on the viscosity to entropy ratio rj/s from hydrodynamical simulations |156] have 
cast strong doubts on the applicability of LO-HTL as a model of the dense medium created at RHIC. 
In addition, a computation of q at NLO in the HTL effective theory was found to introduce a large 
constant term in addition to the log[4£'T/m|)] |157] . This has motivated the authors of Ref. [87] to 
attempt a global comparison to both the light and heavy quark energy loss data assuming that the q, e 
and 62 depend on temperature according to their dimensions, i.e., g oc 62 oc Te oc T^. The longitudinal 
coefficients e and 62 are related by the fluctuation dissipation theorem and they also assumed q = 2e2- 
The one remaining overall constant was set by comparing light quark energy loss to data. With these 
simplifying assumptions, these authors obtained a marked improvement in comparison with the heavy 
flavor suppression. This is illustrated in the bottom panel of Fig. [23] which compares non-photonic 
electron Raa with model curves using a Woods-Saxon overlap geometry with an initial temperature 
of T = 400 MeV at the hottest point. The solid line shows the result using HTL values for the 
transport coefficients q ~ 2GeV^/fm at a jet energy of 20 GeV, e and 62 are calculated using the HTL 
approximation and the same value of as as for q (see also Section [^^2]) . The dashed line shows the result 
without any logarithmic dependence of the transport coefficients on the energy of the jet. The same 
calculation also describes the measured elliptic flow for light hadrons and non-photonic electrons. 

Other models, which include hadron formation inside the medium followed by dissociation and 
recombination |158j or resonant scatterings |159] can also be tuned to describe the measurement. The 
implications of such models for light quark energy loss have not been studied in detail. 

High-pT measurements of heavy flavor production at RHIC have so far relied on electrons from semi- 
leptonic decay and therefore do not distinguish between charm and bottom. The difference in mass 
between charm quarks (m ^ 1.3 GeV/c^) and bottom quarks (m ~ 4.2 GeV/c^) leads to a different 
momentum dependence of the dead-cone effect, which will be explored by future measurements. At 
high pt, the charm quark behaves as a light quark, so that comparisons to light hadrons are a measure 
of the difference between quark and gluon energy loss. Calculations have shown that this effect leads to 
large effects in the heavy-to-light ratios at RHIC and LHC |160] . Future measurements will certainly 
be sensitive to these effects. 

6.4.8 7-hadron and jet based measurements at RHIC 

In recent years, two new types of measurements of parton energy loss are being explored at RHIC: 
7-hadron correlations and jet reconstruction. Both measurements provide experimental control over 
the initial state parton kinematics: in the case of 7-hadron measurements, the photon and the recoil 
jet have equal pt, while in the case of jet reconstruction, the parton energy is measured by summing 
up the energies of the fragments. 

The 7-hadron measurement is conceptually appealing, and has been advocated for a long time 
|165[ 1166] , but it is experimentally difficult because of the low rates and the large background from 
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Figure 24: Recoil hadron suppression for direct photons measured by PHENIX (left panel) |161] and 
STAR (right panel) [162] in central Au+Au collisions compared to p + p. The PHENIX measurement 
(left panel) is presented for four ranges in between 5 and 15 GeV/c and compared to energy loss 
calculations by Zhang, Owens and Wang (ZOWW) |163] . The preliminary STAR measurement (right 
panel) has 8 < < 16 GeV/c and is compared to the ZOWW calculation (dotted curves), a calculation 
by Renk |164] and an AMY-based calculation that takes into account parton-photon conversions in the 
medium (dashed line) [ID] and a similar measurement for 7r° recoil. Figures reproced with permission 
from |161] (©American Physical Society) and |162] . 



7r° decays. First results from PHENIX [161] and STAR (preliminary) |162] are shown in Fig. [2l] and 
compared to model calculations. Both experiments show a recoil suppression of 0.3, which is similar to 
the value found for di-hadron suppression in Fig. [18] (charged hadrons) and the right panel of Fig. IMI 
(yr^-charged hadron pairs) despite the different geometrical and fragmentation bias. 

The 7-hadron measurements are compared to model calculations by three groups. The calculation 
by Zhang, Owens and Wang [163] . which is shown in both panels of Fig. |2l]uses NLO pQCD to calculate 
the initial photon rate, including fragmentation photons but uses a simplified geometry (hard spheres) 
and an energy loss calculation without energy loss fluctuations. The calculation by Renk [164] (solid line 
in the right panel of Fig. (241) uses a more realistic geometry and ASW quenching weights for energy loss. 
The third calculation in the right panel of Fig. [2H (dashed line) uses a 3D hydro model for the medium 
and an AMY-based energy loss formalism, which includes production of photons by interactions with 
the medium |10]. All three model calculations are in good agreement with the measurement, indicating 
that the measured recoil suppression is consistent with the hadron suppression measurements. 

For sufficiently high photon energy, one would expect to see the ^-dependent suppression of the 
fragmentation functions, similar to what is shown in Fig. [TOl although the effect may be weaker for 
a realistic nuclear geometry than for the fixed-length calculations presented in Fig. [TOl Some of the 
model curves in Fig. [24] indeed show a rise at small zt = Pt/Pt ~ 0-2, where no measurements are 
available due to the large combinatorial background at low pt- Future runs at RHIC will bring larger 
data samples and correspondingly reduced uncertainties. 

Experimentally, jet measurements have the advantage of a large cross section, but a good jet analysis 
is complicated, especially in the presence of uncorrelated backgrounds from the underlying heavy ion 
event. Jet measurements are being actively pursued at RHIC, by both STAR [1671 11681 1169] and 
PHENIX [170] . Initially, it was thought that it would be necessary to reject all particles below a certain 
transverse momentum threshold from the analysis to reduce the background, but experience with data 
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has shown that this is not necessary. This enables infrared safe algorithms to be used, such as the 
recently developed SISCone algorithm [171j and fast /ct |172j and anti-Zc^ algorithms |173j . 

In collisions of elementary particles, e"*" + e~, p + p, and p + p, reconstructed jet energies provide 
an accurate measure of the underlying parton momentum as evidenced by kinematic reconstruction 
of W ^ qq \174:\ 11751 1176] and t — >■ bW qq decays \177\ I178j . The accurate mapping between 
parton energy and jet energy is a phenomenological implication of QCD factorization: an infrared 
and collinear safe jet definition is insensitive to details of the soft QCD processes in jet fragmentation 
and therefore measures properties of the 'hard partons'. In heavy-ion collisions, jet-fragmentation is 
modified by interactions with the medium, thereby modifying the longitudinal and transverse structure 
of the jet. In addition, part of the jet energy may be transported outside the jet cone by large angle 
radiation, thus reducing the reconstructed jet energy. A number of different experimental analyses are 
being performed, which are sensitive to different aspects of the medium modification of jet energies. 
The available preliminary results show an increasing suppression when using a smaller jet radius in the 
reconstruction algorithm, suggesting that jet broadening is an important effect |168j . For large radius, 
R = 0.4, the jet suppression is much smaller than the inclusive hadron suppression, indicating that a 
large fraction of the radiated energy is recovered by jet reconstruction. 

Attempts have also been made to measure fragmentation functions, i. e. the longitudinal distribution 
of particles in jets, with the aim to characterize the softening of the fragmentation due to medium in- 
teractions. Such measurements, however, are complicated by potential interplays between the softening 
of the fragmentation and the jet energy reconstruction. Preliminary results that attempt to circumvent 
this problem by using di-jet events, where the jet energy is measured in one jet and the fragment dis- 
tribution in is studied in the recoil jet, show a large suppression of the overall recoil jet rate, but no 
apparent modification of the fragment distribution in the reconstructed recoil jets |179j . 



6.5 Outlook: how to reduce uncertainties 

In the above various measurements have been compared to theoretical calculations to extract medium 
parameters. In some cases, such as the heavy flavor measurements, the experimental uncertainties are 
sizable, of the order of 20% and seem to limit the extraction of medium properties. As mentioned 
in Section 16. the first approach to improving this situation should be by defining more selective 
observables using theoretical insight and modeling. For example, one might expect a measurement of 
recoil suppression Iaa of heavy quarks to be sensitive to the path length dependence of energy loss, like 
in the case of light hadrons. This is clearly a challenging measurement, so it would be worth to first 
explore the potential in model calculations. 

However, it is also worth to consider the origins of the various experimental uncertainties to see 
whether improvements are possible. In this section, we will briefly discuss the various sources of uncer- 
tainties in the measurements and which improvements are planned with future RHIC upgrades. 

Statistical errors, luminosity 

The most straightforward source of uncertainties are statistical uncertainties. These uncertainties 
are always due to counting statistics in the experiment and can therefore be reduced by increasing the 
number of analyzed events, either by improving luminosity or by longer running times, and sometimes 
by selecting a process with a larger cross section. The beam luminosity and effective running time at 
RHIC have steadily increased over the years. The largest event samples were taking in run 4 and run 7, 
with run 7 having about twice the integrated luminosity of run 4. It is important to realize that not all 
observables can make full use of the luminosity, because the event rates may be too high to record all 
events. Both STAR and PHENIX can select events with large energy depositions in the electromagnetic 
calorimeters with an on-line trigger, so that measurements with high-p^ photons, 7r° and electrons can 
make use of the full luminosity. Both experiments are also continually improving the detector read-out 
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rates to be able to record a larger fraction of the total luminosity for other (un-triggered) measurements. 

If a signal has a large background, the statistical uncertainty is driven by the amount of background. 
This explains the relatively large statistical uncertainties on direct photon measurements (decay pho- 
ton background) and to a lesser extent the non-photonic electron measurements (conversion electron 
background) . 

Background uncertainties 

Apart from the statistical uncertainties on a background, there are also systematic uncertainties. 
For example, in the case of the decay photon background for direct photons, the background level is 
set by the 7r° spectrum, which is known to finite precision. The contributions of higher-mass mesons 
are less precisely known. The main source of the conversion electron background for the (non-photonic) 
electron spectrum are also photons from 7r° decay, so here again the subtraction relies on the knowledge 
of the TT*^ spectrum, in addition to a precise knowledge of the amount and composition of the detector 
materials. 

If the signal is smaller than the background, a relatively small uncertainty in the background level 
may translate into a large uncertainty on the signal. Larger data samples can improve the knowledge 
of the background level, but only up to the level of the intrinsic systematic uncertainties. 

Intrinsic systematic uncertainties 

Every experimental set-up has intrinsic lower limits to the experimental uncertainties. The two 
dominant effects are uncertainties in the detection efficiency and in the energy or momentum scale. 
Experimental efficiency is generally calculated using detector response simulations, which have a finite 
precision. Details of what limits the uncertainty differ from detector to detector, but it is safe to 
say that a precision of better than a few per cent in the detection efficiency is difficult to attain. In 
addition, energy scale calibration uncertainties may introduce significant uncertainties in the yields at 
a given p^. For example, with a power law spectrum an energy scale uncertainty of 0.5% leads 

to a 3% uncertainty on the yield. These types of uncertainties limit, for example, the tt^ spectrum at 
intermediate pt in Figs [I6] and [T71 

Normalization uncertainties 

For spectra and Raa, normalization uncertainties are a significant contribution, see Figs [121 [161 ^i^nd 
[23l where this type of uncertainty is indicated by a separate error band without data point. There 
are two main contributions: the absolute normalization of the p+p reference measurement and the 
determination of the number of binary collisions Ncoii- 

The normalization uncertainty of p+p reference measurements is due to a significant contribution 
of low-multiplicity events, including double-diffractive collisions, to the total cross section. Some of 
these low-multiplicity events are not detected by the experiment. The measurement in p+p can be 
given in terms of multiplicity per event {(PN / dydpT) or differential cross section (Pa / dydpT- The 
normalization uncertainty is slightly different for these two types of measurement; in one case it enters 
as an uncertainty in the total number of events, while in the other case it enters as an uncertainty in 
the total cross section, which can be be measured independently, using a Vermeer scan. Both methods 
lead to a total uncertainty on the absolute cross section of 5-10%. 

To calculate Raa oi Rcp in heavy ion events, the number of binary collisions (or overlap function 
Taa) is used. This quantity is derived from Glauber models of the total cross section (cf Section ISTTl) . 
The uncertainty on this calculation is indicated separately in Fig. [16] as a grey band on the left side of 
each panel. The result is accurate for central events, but for peripheral collisions the uncertainties can 
be sizable. This uncertainty can possibly be reduced with dedicated studies of peripheral events, but 
some of the uncertainty is inherent in the Glauber modeling. 

The normalization uncertainty is strongly reduced when concentrating on central events, where the 
determination of N^oii is quite reliable. The types of normalization uncertainty discussed here only affect 
measurements that involve an absolute cross-section and therefore do not affect di-hadron measurements 
and elliptic flow. 
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When comparing theory calculations to experimental data, all of the above sources of uncertainties 
need to be taken into account. For example, in Fig. [T51 there are three types of uncertainties indicated, 
statistical, systematic and an overall normalization uncertainty. The systematic uncertainties are some- 
what correlated between points and are due to the sources listed above under background uncertainties 
and intrinsic uncertainties. For the vr" spectrum in the left panel of Fig. [12], statistical uncertainties are 
dominant only for px > 16 GeV/c. At lower pt, the systematic uncertainty of about 5-10% on each data 
point and the overall uncertainty of about 10% are dominant. Reducing the systematic uncertainties is 
very difficult, and often requires building new detectors. 

For the non-photonic electron measurements in the right panel of Fig. [12] and in Fig. |23l the 
statistical uncertainties are dominant for p^- > 5 GeV/c, so that longer running times or larger luminosity 
will still improve the precision of the result. Note also that the statistical uncertainties for the STAR 
measurement are larger, due to a smaller analyzed (triggered) data sample. Also for this measurement, 
reducing the systematic uncertainties below the current value of 15 — 20% is difficult to impossible. The 
precision of the di-hadron suppression and the 7-hadron measurements are still limited by statistical 
uncertainties, so that future runs with larger luminosity will improve the precision. 
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7 Outlook: parton energy loss at LHC 



In the near future, the LHC will start colliding protons and nuclei at energies far greater than RHIC 
energies (maximum energies are ^/sNN = 5.5 TeV for Pb+Pb and ^/s = 14 TeV for p+p, but the initial 
runs will be at lower energy). 

First results from the LHC experiments should already provide an important test of the understand- 
ing of parton energy loss developed at RHIC: if the same models and theory can be applied to describe 
the LHC results, that would be an important sign that our basic understanding of parton energy loss 
is correct. However, the LHC will provide access to a far greater number of jet observables. 

Our current understanding of RHIC results indicate that jet energy loss is typically a few GeV with 
important tails to much larger energy losses. Since typical jet energies at RHIC are 10-20 GeV, the 
energy loss is not much smaller than the jet energy, which complicates the theoretical treatment of 
energy loss. At LHC, large numbers of jets with energies above 100 GeV will be produced, while the 
medium density is only expected to increase by a modest factor of 3-5 |180] . This makes it possible 
to cleanly separate jets from the underlying event. In addition, the regime where AE < E, where 
theory calculations are better controlled, may become experimentally accessible. Luminosity upgrades 
at RHIC will also extend the accessible pr-Tange there, reaching parton energies of 40-50 GeV. 

The main jet observables have already been discussed in the context of RHIC results in Section 
16.4.81 The basic measurement is the jet cross section which shows to what extent the radiated energy 
is recovered by jet-finding algorithms. Once this is established, measurements of longitudinal and 
transverse jet structure (fragmentation function and energy profile ^'(-R)) are used to quantify the 
redistribution of energy inside the jet. An example calculation of modified fragmentation functions in a 
fixed-length medium is given in Fig. [TOl However, jet reconstruction intrinsically involves multi-hadron 
observables, which are best modeled using Monte Carlo event generators, such as discussed in the next 
section. The large jet rates at 100 GeV at LHC will further enable more differential studies of e.g. jet 
structure using sub-jet distributions [181] . 

7.1 Experimental capabilities at the LHC 

ALICE |182[ 1183] is the dedicated heavy-ion detector, which has been optimized for large track densities 
in heavy-ion collisions and features powerful particle identification. ATLAS [184[ I185j and CMS |186] 
will also perform heavy-ion measurements with their detector systems that are more conventional for 
high-pj- and jet measurements, including large calorimeter coverage. For a complete review of the 
potential of hard probes at the LHC, see |124[ 1187] . 

The LHC detectors are all equipped with high-precision vertex detectors for secondary vertex re- 
construction of heavy flavor decays. This will allow to separate charm and beauty from the outset, 
providing critical tests of the mass and flavour dependence of energy loss, which will show whether 
radiative energy loss is dominant or not. Since the dead-cone effect should be less pronounced for the 
lighter charm quark than for the heavy bottom quark, comparing the suppression of charm and bottom 
provides direct sensitivity to the dead-cone effect. On the other hand, comparing charm fragmentation 
at high pt with light parton fragmentation, which is a mix of quark and gluon fragmentation, provides 
sensitivity to the difference between quarks and gluons |160] . In addition, heavy quark fragmentation 
studies should be possible at LHC, but these have not yet been explored in detail. 

Finally, 7-jet and Z^-^et events will be studied to provide a well-defined calibration of the jet 
energy measurement. The total rates of these events are significantly lower than inclusive jets, so 
that differential measurements of jet fragmentation will most likely be performed with single jets or 
di-jet events. For a more detailed discussion of the potential of such measurements, and the expected 
performance of the LHC detectors, see |188[ 1186] . 
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7.2 Theoretical issues: From analytical approaches to Monte-Carlo 

While the single hadron formalism has been the central theme in this review we now report on the 
emerging new directions in jet modification theory which are being tailored to calculate multi-hadron 
final states. 

The computation of multihadron final states for jets in vacuum is now a rather well established 
procedure. The standard method uses Monte-Carlo routine to compute the underling jet shower by 
sampling the DGLAP splitting functions repeatedly in order to sequentially decompose a single virtual 
parton into a shower of lower virtuality partons, using the Sudakov prescription. We refer the reader 
to Refs. |189[ 1190] for details on these algorithms. At any given step in this routine one typically has 
a hard virtual parton which will survive for a time which is inversely proportional to its virtuality and 
then split into two partons with lower virtuality and so on. At low enough virtuality, this "DGLAP" 
shower routine is terminated and a phenomenological hadronization routine is invoked to convert the 
parton shower into a shower of hadrons. 

In a medium, the showering and hadronisation process becomes far more complicated, due to addi- 
tional radiations and the possible exchange of color with the medium |191] . In addition, hadronisation 
through recombination of shower partons with partons from the medium might play a role. At the 
LHC, given the energies of the jets, one expects a large part of the jet to exit the medium in partonic 
form and continue to shower. Finally one expects a considerable part of this escaped and modified jet 
to hadronize independently of the medium. While the very soft part of this hadronic spectrum from 
the jet will still become indistinguishable from the harder hadrons from the medium, the high px part 
of the jet should be separable from the medium. The major challenge is to model the partonic split- 
tings which now contain a vacuum like part where a higher virtuality parton splits into lower virtuality 
partons interfering with an in-medium part where scatterings in the medium may raise the virtuality 
of the hard parton and cause it to split into two. 

One of the first Monte-Carlo event generators to include an implementation of parton energy loss 
was PYQUEN/HYDJET [192] . The model includes elastic and radiative energy loss. The radiative 
energy loss is based on the BDMPS formalism and the radiated gluons are added to the partonic final 
state, which is then fragmented using PYTHIA 6.2. HYDJET also generates a soft underlying event 
using a simple parametrisation including elliptic fiow. The programs is also available in a C-I--I- version 
HYDJET++ [193], which is still being developed. 

Besides PYQUEN, there are currently four different programs: JEWEL [HI], Q-PYTHIA [M], 
YaJEM |195] . and MARTINI |196] in the published literature, which have been benchmarked with single 
hadron data. Most of these generators are also still under active development by the authors. They all 
use a somewhat different set of approximations and calculational techniques. Of these Q-PYTHIA |194 
is closest in spirit to the traditional vacuum Monte-Carlo routines. Here the authors add a medium- 
induced contribution to the vacuum splitting function and use this to calculate an in-medium Sudakov 
form factor. They now have a probability that the jet will emit a gluon while propagating over a certain 
distance and also the probability that it will not radiate. This is then iteratively applied to compute a 
shower distribution. In the case of JEWEL |181] . the authors attempt an original algorithm. The very 
process of single gluon radiation is calculated via a numerical algorithm. Here, the random process of 
scattering is used to modify the formation time of the radiated gluon in a stochastic process. More 
scatterings tend to increase the transverse momentum of a gluon which has not completely formed. As 
the transverse momentum is increased the formation time of the gluon shrinks till it is formed within the 
next time step. After this the process begins anew. JEWEL also includes the elastic loss that the hard 
jet encounters every time it scatters. In YaJEM |195] . the author modifies the vacuum DGLAP shower 
that is contained within PYTHIA. This is done by increasing the virtuality of an intermediate state 
due to scattering in the medium. This leads to a change in the radiation pattern and thus to a different 
shower pattern. In MARTINI |196] . analytical energy loss routines computed within the AMY-McGill 
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formalism including elastic energy loss are used to modify the shower distribution in PYTHIA prior to 
hadronization. MARTINI uses tabulated hydrodynamic data for the medium density evolution, with a 
choice of hydro dynamical results from different groups. 

These Monte-Carlo routines along with others being developed represent the next wave of jet mod- 
ification calculations. These will no longer be limited to the single hadron or dihadron level but will 
be able to compute a whole host of full jet observables. They will also provide the experimental com- 
munity with multiple event generator codes which will provide necessary input and benchmarks for the 
experimental analysis of multihadron final states. 
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8 Conclusion 



In this review, we have discussed in some detail the theoretical formalisms describing parton energy 
loss in cold and hot QCD matter. Conceptually, there are two main mechanisms for energy loss: elastic 
and radiative. For radiative energy loss off high-momentum partons in a dense medium, the formation 
times are large, so that the Landau-Pemeranchuk-Migdal (LPM) effect plays an important role, leading 
to energy loss scaling with the square of the path length (L^). At the same time, a high-momentum 
parton is also expected to reduce its virtuality even in the absence of a medium, leading to a jet of 
hadrons. The interplay (and interference) between vacuum radiation and medium-induced radiation is 
therefore also expected to play an important role. 

Due to the complex nature of the energy loss process and in-medium modification of the fragmenta- 
tion function, different theoretical formalisms have chosen different approximations in the calculation. 
A side-by-side comparison of the medium modifications to the fragmentation functions for a fixed-length 
homogeneous medium was presented in Fig. [10] and shows important quantitative differences between 
different formalisms. 

Unfortunately, fixed length, homogeneous hot QCD matter is not experimentally available, but 
instead a rapidly expanding medium is produced in heavy ion collisions. Comparing Fig. [T7] and [TUl we 
concluded that the way in which the partons sample the evolving medium is an important aspect of the 
modeling for heavy ion collisions, which currently precludes the precise determination of the medium 
density or transport coefficient from experiment. 

There are, however, qualitative features of radiative energy loss that can be addressed with ex- 
perimental data despite the uncertainties in some of the quantitative details of the calculations. The 
comparison of measured nuclear modification factor Raa and the recoil suppression I^a indicate that 
there is a leading dependence to the energy loss (cf. Fig. 12^ . which is a first indication that the 
LPM interference effect indeed plays an important role in parton energy loss. 

Another prediction of radiative energy loss is that heavy quarks lose less energy than light quarks due 
to the dead-cone effect. The measured suppression of heavy quarks at RHIC is larger than expected from 
pure radiative energy loss, indicating that other effects play a role. Possibilities include the presence of 
elastic energy loss and the possible influence of the larger scale at which transport coefficients should 
be evaluated due to the larger mass of the heavy quark. 

Future measurements at RHIC and LHC will include 7-jet measurements and measurements with 
kinematically reconstructed jets, which will provide more differential control of the partonic kinematics. 
Such measurements will be more directly comparable to the theoretically calculated fragmentation 
functions in Fig. [TUJ 

An important ongoing theoretical development is the use of Monte Carlo techniques, which has two 
important aspects. On the one hand, Monte Carlo techniques can be used to explore the effect of certain 
limitations of the analytical calculations of energy loss, most notably the effect of kinematic limits. On 
the other hand, Monte Carlo event generators can be used to calculate more exclusive (multi-hadron) 
observables, and thus provide a reliable tool to explore observables based on reconstructed jets, such as 
fragmentation function and jet broadening measurements. 

All in all, a quantitative understanding of the theory and phenomenology of parton energy loss 
based on the RHIC measurements is starting to emerge, including some areas of apparent disagreement 
between the data and theory. As a result, all the tools and understanding are available to start devis- 
ing critical tests of the theoretical understanding using future measurements at RHIC and LHC. The 
much larger momentum range available at LHC will likely make the regime where the energy loss is 
significantly smaller than the jet energy experimentally available. This regime is the most suited for 
pQCD based theoretical treatment and therefore promises to provide the most unambiguous tests of 
our understanding of fragmentation and energy loss in a hot QCD medium. 
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